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all-electron  beam  storage  concepts  (alloy  junction  and  electron- 
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PREFACE 


The  objective  of  this  basic  research  program  is  to  demonstrate 
the  feasibility  of  an  Advanced  Archival  Memory  System  capable  of 
permanent  storage  of  data  on  solid  planar  media  with  0.1  micron 
bit  spacing  and  10  MHz  read  and  write  rates.  The  eventual  goal 
of  this  program  is  the  construction  of  low-cost  memory  systems 
of  lO^^-lO^®  bit  size  with  megabit  per  second  data  rates  and  fast 
access  times.  The  data  will  be  stored  on  semiconductor  target 
plates  by  milling  holes  (bit-site)  in  a surface  layer  of  a selenium- 
arsenic  mixture  using  a high  current  focused  electron  probe.  The 
bit  pattern  is  read  by  passing  a low  current  focused  electron  probe 
over  the  target.  Vfhere  there  is  a milled  hole,  the  electron  beam 
will  penetrate  into  the  silicon,  and  the  resultant  current  will 
be  amplified  by  a diode  field  structure  integral  to  the  target. 

Where  there  is  no  hole,  no  signal  will  be  picked  up.  The  modula- 
tion of  the  current  through  this  diode  as  the  electron  beam  moves 
from  one  bit  site  to  another  provides  the  read  signal. 
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This  report  covers  the  effort  for  the  first  eight  months  of 
1978.  It  is  the  second  report  in  Phase  II  of  this  program.  Phase  I 
was  covered  in  Reports  AFAL-TR-76-213,  AFAL-TR-77-35,  AFAL-TR-77-157, 
and  AFAL-TR-77-202.  In  this  reporting  period,  the  contract  has  been 
redirected  from  an  ion  damage  writing  concept  to  the  electron  beam 
milling  approach  described  above.  Background  surrounding  this  re- 
direction can  be  found  in  the  introduction  to  this  report  and  in 
the  first  Phase  II  report,  AFAL-TR-78-96.  The  work  during  this 
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period  has  concentrated  on  evaluating  two  different  all-electron 
beam  storage  concepts  (alloy  junction  and  beam  milling)  and  choos- 
ing a preferred  approach  for  further  work.  The  work  has  shown 
that  the  alloy  junction  approach  is  not  presently  feasible  and 
that  the  electron  beam  milling  approach  looks  quite  feasible.  In 
place  of  the  ion  beam  writer  that  was  to  be  built,  a high  current 
electron  probe  is  being  purchased  for  use  in  evaluating  the  elec- 
tron milling  storage  concept. 
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SECTION  I 
INTRODUCTION 


BACKGROUND 

The  objective  of  this  basic  research  program  is  to  investi- 

14  15 

gate  the  feasibility  of  new  storage  techniques  for  10  - 10 

bit  archival  memories  using  ion  and  electron  beams  with  semicon- 
ductor targets.  The  current  effort  (Phase  II)  deals  with  1)  the 
demonstration  of  the  feasibility  of  writing  and  reading  at  data 
rates  of  10  megabits  per  second  with  a storage  density  of  0.1  ym- 
bit  size,  2)  the  development  of  the  storage  medium,  and  3)  the 
initiation  of  system  application  studies. 

Three  candidate  storage  methods  have  been  or  are  being  evalu- 
ated. The  first  of  these  is  the  formation  of  small  damaged  re- 
gions on  the  surface  of  silicon  diode  detectors  with  a focused 
ion  beam.  The  other  two  approaches  are  based  on  the  use  of  elec- 
tron beam  heating  of  metal-silicon  target  structures.  The  first  ! 

of  these  methods  is  based  on  forming  surface  diodes  by  heating 
a thin  metallic  coating  on  a silicon  target  above  the  silicon-  I 

metal  eutectic  temperature  with  a finely  focused  electron  bean.  j 

The  goal  is  to  regrow  sufficient  silicon  as  the  mixture  cools  | 

down  such  that  the  silicon  will  then  be  doped  by  the  metal  specie  I 

and  a microscopic  diode  is  formed.  The  second  approach  is  either 
to  vaporize  or  to  deform  an  amorphous  metal  layer  with  the  electron 


bean  so  that  the  written  area  will  transmit  an  electron  beam  when 
read,  while  the  unwritten  area  will  not.  Conceptually  this  is 
a very  simple  idea.  The  metal  layer  is  essentially  a punched  card 
that  is  read  by  an  electron  beam.  Unlike  the  previous  two  approaches, 
the  silicon  is  not  altered  in  the  writing  process  and  only  func- 
tions as  a detector  for  the  electron  beam.  The  goal  for  all  three 
approaches  is  to  store  and  retrieve  bits  at  a 0.1  micrometer-density 
at  a rate  of  10  megabits  per  second. 

The  program  has  been  redirected  during  this  period  of  the 
contract.  Originally  the  contract  called  for  building  a focused 
ion  writing  optical  system  and  an  electron  beam  reading  system. 
However,  these  items  will  not  be  carried  out.  In  their  place  a 
state-of-the-art  high  current  electron  writing  column  will  be 
purchased  for  use  in  evaluating  the  preferred  electron  beam  writing 
approach.  The  change  in  direction  in  the  contract  came  as  a re- 
sult of  determining  that  the  ion  writing  approach  would  fail  in 
the  stated  goal  of  O.lwm-bit  density.  The  following  is  a brief 
description  of  the  conclusions  in  this  area.  A complete  discus- 
sion is  found  in  the  Phase  II  Interim  Report  (July  15-December 
31,  1977). 


In  using  writing  fluences  consistent  with  the  write  column 
design,  the  resolution  appeared  to  be  limited  to  0.30  um.  The 
bit  resolution  criterion  used  was  equal  signals  from  isolated 
written  and  unwritten  areas.  This  is  a stringent  but  correct 
criterion.  Features  close  to  0.1  pm  in  size  could  be  easily  iden- 
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tified  from  oscillograph  readouts,  but  simple  feature  definition 
is  not  adequate  basis  for  a memory  technology.  Considerable  effort 
was  also  spent  in  determining  the  factors  limiting  resolution. 

Two  resolution  limiting  effects  were  found.  These  include  a defect 
diffusion  phenomenon,  in  which  written  defects  migrate  beyond  their 
original  sites  and  a beam  focusing  problem.  In  the  case  of  the 
latter,  the  focused  electron  beam  used  in  the  experiments  has  a 
highly  concentrated  component  of  electron  current  in  the  beam 
center  which  determines  the  ultimate  resolution  for  secondary  elec- 
tron beam  microscopy.  It  also  appears  to  have  had  wide  tails 
which  contained  sufficient  current  to  produce  much  poorer  overall 
resolution  for  an  archival  memory. 

The  work  in  the  current  reporting  period  has  concentrated 
on  evaluating  two  approaches  to  an  all  electron  beam  based  memory, 
alloy  junction  and  electron  machining.  Two  additional  tasks  are 
also  being  carried  out.  The  first  is  a read  systems  study  to  an- 
swer questions  regarding  the  size  of  the  bit  field  that  can  be 
electronically  accessed  and  the  compatibility  of  reading  and  writ- 
ing with  the  same  column.  The  second  is  an  ion  source  test  that 
will  complete  the  design  and  specifications  of  the  ion  column  be- 
gun earlier  in  Phase  II.  In  the  case  of  the  latter,  the  ion  sources 
obtained  from  the  Oregon  Graduate  Center  are  being  tested  for  angu- 
lar brightness  and  energy  spread  using  a test  stand  now  completed. 

Both  approaches  to  the  memory  technology  use  below  board 
silicon  diode  detection  for  readout.  The  alloy  junction  work 


is  a continuation  of  the  work  in  this  area  started  in  Phase  I. 

The  electron  machining  approach  builds  on  work  done  in  the  Elec- 
tronics Laboratory  of  the  General  Electric  Company  in  1968-1971 
and  has  made  maximum  use  of  the  results  achieved.  The  efforts 
at  the  Electronics  Laboratory  failed  due  to  problems  with  target 
fabrication,  with  the  secondary  readout  used,  and  with  controlling 
and  getting  sufficient  current  from  the  electron  writing  column. 

To  achieve  sufficient  readout,  the  amorphous  metal  layer  was  made 
from  a triple  alloy  system  of  selenium,  arsenic,  and  gold.  It 
was  difficult  to  control  the  composition  of  these  targets,  and  the 
resolution  achieved  was  quite  variable.  In  addition,  the  beam 
power  available  at  the  time  was  not  sufficient  to  vaporize  the 
selenium  and  the  writing  was  based  on  melting  the  selenium  and 
leaving  a depression. 

Now  over  six  years  of  work  on  diode  detection  in  the  BEAMOS 
project  and  Phase  I and  II  of  this  contract  at  General  Electric 
Corporate  Research  and  Development  offer  promise  of  eliminating 
the  reading  difficulties.  In  addition,  a great  deal  of  progress 
has  been  made  in  both  controlling  field  emission  electron  sources 
and  in  obtaining  higher  overall  currents.  An  order  of  magnitude 
of  improvement  has  been  made  in  the  latter  since  1969.  As  a re- 
sult, it  is  now  possible  to  write  by  vaporizing  the  film.  High 
signal-to-noise  ratios  are  achieved  because  of  the  effective  ma- 
chining and  the  inherent  amplification  occurring  with  diode  read- 
out. The  target  manufacture  is  basically  simple  because  only  a 
selenium  arsenic  mixture  is  used  in  the  film  and  the  film  proper- 
ties are  not  sensitive  to  the  detailed  composition. 

4 


1^ 


CURRENT  STATUS 


An  initial  series  of  exper iements  on  the  electron  beam  milling 
of  selenium  arsenic  target  structures  has  been  performed  and  ana- 
lyzed. A theory  has  been  developed  which  correctly  explains  the 
experimental  results  and  has  been  used  to  predict  the  necessary 
conditions  of  beam  current,  beam  voltage,  spot  size,  and  film  thick- 
ness to  achieve  0.1  pm-bit  densities  and  a 10  MHz  writing  rate. 
These  conditions  have  been  shown  to  be  consistent  with  the  spec- 
ifications for  an  advanced  field  emission  electron  column  being 
purchased  from  ETEC  Corporation.  The  decision  to  purchase  this 
column  was  made  after  a careful  build  or  buy  studying  process. 
Ultimately  the  decision  was  made  because  the  predicted  performance 
of  the  ETEC  microscope  comes  very  close  to  the  optimum  achievable. 

Work  on  the  alloy  junction  writing  approach  has  also  progressed 
a great  deal  in  this  reporting  period.  A detailed  steady  state 
model  of  the  alloying  junction  writing  process  has  been  developed. 
The  model  can  be  used  to  choose  the  alloy  system  which  will  pro- 
duce the  largest  writing  efficiency  and  also  to  determine  the  con- 
ditions of  temperature  and  metal  thickness  necessary  to  achieve 
the  writing.  The  model  was  then  applied  to  the  specific  cases 
of  the  gold-silicon  and  the  aluminum-silicon  alloy  systems.  In 
addition,  a series  of  steady  state  alloying  experiments  have  been 
carried  out.  These  experiments  had  the  desired  result  of  confirm- 
ing the  correctness  of  the  alloying  model  but  also  clearly  iden- 
tified problems  with  the  alloy  junction  approach.  One  significant 


problem  is  that  temperature  excursions  of  over  300  have  been 
theoretically  and  experimentally  shown  to  be  necessary  for  adequate 
writing,  but  the  modeling  of  the  electron  beam  heating  effects 
show  that  excursions  will  be  limited  to  under  75  for  both  pre- 
ferred alloy  systems  (Al-Si  and  Au-Si) . These  temperature  excur- 
sions were  calculated  on  the  basis  of  0.1  mA  current  in  a 0.1 ym- 
spot.  This  is  the  upper  end  of  the  expected  performance  from  the 
state-of-the-art  electron  column  being  purchased  from  ETEC  Corpora- 
tion and  will  not  be  materially  improved  upon  in  the  foreseeable 
future.  In  addition  to  the  difficulty  in  reaching  the  required 
temperature  rises  a problem  of  non-wetting  of  the  metal  specie 
on  the  silicon  surface  was  discovered.  When  heated,  the  metal 
specie  has  a tendency  to  move  away  from  the  heat  source  rather 
than  form  an  alloy.  There  is  also  a tendency  for  the  metal  to 
condense  on  surface  faults  on  the  silicon  target. 

The  electron  milling  approach  has  been  chosen  as  the  preferred 
approach.  This  decision  was  arrived  at  because  of  the  difficulties 
with  the  alloy  junction,  on  one  hand,  and  the  successful  results 
and  optimistic  predictions  for  the  electron  milling  approach,  on 
the  other  hand.  Work  will  continue  on  furthering  our  understand- 
ing of  the  writing  mechanism  and  in  experimenting  with  and  obtain- 
ing good  quality  selenium-arsenic  films.  Upon  arrival  of  the  ETEC 
column,  a program  will  be  started  to  demonstrate  the  storage  capa- 
bilities of  electron  milling  of  selenium  films.  This  will  include 
achievable  bit  density  and  writing  rate. 
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Construction  of  the  ion  source  test  stand  has  been  completed 
and  it  is  now  fully  operational.  Energy  spread  and  angular  bright- 
ness data  are  being  taken  on  the  field  ion  sources  from  the  Oregon 
Graduate  Center.  Results  of  these  measurements  will  be  included 
in  the  final  report  under  the  current  contract.  Work  on  the  read 
systems  study  to  investigate  the  limits  on  the  size  of  directly 
addressable  data  blocks  and  on  the  use  of  matrix  lens  technology 
will  be  started  in  the  last  quarter  of  calendar  1978. 

REPORT  CONTENTS 


In  Section  II  of  the  report  the  steady  state  alloying  model 
is  presented  in  detail.  Specific  examples  of  the  Al-Si  and  Au— Si 
alloy  systems  are  presented.  Experimental  data  on  the  corrobora- 
tion of  the  model  using  an  Au-Si  alloy  junction  are  also  presented. 
Section  III  discusses  the  results  of  isothermal  studies  of  the 
alloying  process.  The  problems  of  wetting  the  silicon  surface 
by  the  metal  specie  during  alloying  are  demonstrated  and  discussed. 

In  Section  IV  of  this  report  the  results  of  the  selenium- 
arsenic  writing  experiments  and  the  theoretical  analysis  are  pre- 
sented in  detail.  A computer  program  has  been  developed  to  calcu- 
late transient  heating  of  target  structures  in  two  dimensions  with 
up  to  three  distinct  layers.  This  model  is  used  in  the  electron 
milling  and  the  alloy  junction  analysis  and  is  fully  explained 
in  Section  V of  the  report.  The  program  and  a user's  manual  are 
presented  in  the  Appendix.  Analysis  of  the  ETEC  microscope  is 


I presented  in  Section  VI  where  the  predicted  performance  of  the 

electron  optical  column  is  shown  to  agree  well  with  the  published 
i:  specifications  and  data.  The  phenomena  limiting  performance  are 

f 

also  discussed. 


SECTION  II 

SILICON-METAL  ALLOY  MODEL 


INTRODUCTION 


During  this  reporting  period  consideration  was  given  to  the  ^ 

alloy  process  for  binary  phase  systems.  Based  on  a general  binary  | 

system  phase  diagram,  a model  was  developed  that  explains  the  alloy 
process  from  a steady  state  point  of  view.  Use  of  the  model  has 
aided  in  understanding  the  requirements  and  conditions  for  a good 
binary  archival  target  system.  The  model  has  also  been  useful 
for  choosing  laser  writing  experiments.  The  development  and  use 
of  the  alloy  model  during  this  reporting  period  are  discussed  in 
this  section. 

THE  ALLOY  MODEL 


1 i 

t 

! 
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Many  of  the  metals  and  dopants  found  in  silicon  semiconductor 
processing  exhibit  a binary  phase  with  silicon.  The  phase  dia- 
gram for  these  materials  is  of  the  type  shown  in  either  Figure  la 
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or  Figure  lb.  In  Figure  la  the  system  shows  a eutectic  tempera- 
ture (Tg)  that  is  depressed  from  the  melting  points  of  either  con- 
stituent and  a reasonable  weight  percentage  (Wg)  of  silicon  in 
the  eutectic  melt.  Typical  materials  showing  this  type  of  phase 
diagram  are  Al,  Au,  and  Ag;  and  are  representative  of  metals  with 
moderate  melting  temperatures.  Figure  lb  shows  a system  where 
the  melting  point  of  the  material  interacting  with  silicon  and 
the  eutectic  temperature  are  almost  the  same.  Of  further  note 
for  these  systems  the  weight  percent  of  silicon  in  the  eutectic 
melt  is  extremely  small.  Materials  exhibiting  this  behavior  are 
typically  of  the  low  melting  point  metals  such  as  Ga  and  In. 

The  alloy  process  can  be  explained  by  referring  to  the  phase 
diagrams  of  Figures  la  and  lb  and  the  material  cross  section  of 
Figure  2.  This  figure  shows  a silicon  substrate  with  a thin  metal 
layer  and  a small  alloy  region  presumably  formed  by  heating  the 
material  at  that  point  with  localized  radiation  such  as  an  elec- 
tron or  laser  beam.  For  purpose  of  explanation,  a metal  over-layer 
is  assumed  because  binary  phase  diagrams  seem  to  be  more  typical 
of  metal-silicon  systems.  However,  it  is  acknowledged  that  metal- 
loids (e.g.  Sb)  can  have  binary  phase  with  silicon  and  other  con- 
ventional dopants  such  as  As,  P,  and  B show  a binary  type  phase 
over  a portion  of  their  phase  diagrams.  The  following  discussion 
and  the  resulting  model  apply  equally  well  to  these  systems. 

Referring  now  to  the  alloy  region  of  Figure  2,  the  region 
was  formed  by  raising  the  local  temperature  to  a maximum  value 
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Figure  1.  Typical  Binary  Phase  Diagrams  for  Si 
and  materials  Used  in  Si  Processing 


which  is  greater  than  the  eutectic  temperature  as  indicated 
on  Figure  1.  It  is  assumed  that  this  temperature  is  maintained 
and  relaxed  at  rates  so  that  equilibrium  conditions  are  met.  The 
resulting  alloy  structure  that  is  formed  consists  of  a regrown 
crystal  region  of  thickness,  tj^,  and  a homogeneous  eutectic  region 


METAL 


of  thickness,  t^.  The  reason  that  the  regrown  crystal  forms  is 
that  the  solubility  of  silicon  in  the  alloy  melt  is  reduced  as 
the  temperature  is  lowered  from  Tj^  (see  Figure  1)  . With  a finite 
thickness  of  metal  present,  the  melt  region  must  reduce  in  size 
as  the  temperature  is  lowered,  thus  leaving  a regrown  crystal. 
These  statements  will  be  quantified  below.  The  regrown  crystal 
region  is  doped  with  metal  according  to  the  solid  solubility  of 
the  metal  in  silicon  at  the  temperature  at  which  the  region  passes 
through  the  solidus.  Solid  solubilities  [1]  as  a function  of 
temperature  of  some  of  the  typical  dopants  used  in  silicon  sem- 
iconductor processing  are  shown  in  Figure  3.  As  indicated  on  the 


phase  diagram  of  Figure  1 the  alloying  process  follows  the  phase 
diagram  on  the  silicon  rich  side  of  the  eutectic  since  there  is 
an  overabundant  supply  of  silicon  in  comparison  to  the  thin  metal 
layer . 

The  depths  to  which  the  alloy  regions  reach  are  therefore 
determined  by  the  temperature  and  weight  percentages  of  the  phase 
diagram  as  well  as  the  metal  layer  thickness.  The  alloy  depths 
can  be  quantified  by  considering  the  weight  percentage  (Wj^j)  of 
silicon  in  the  melt  formed  by  heating  to  a maximum  temperature 


where 


W, 


M 


Ps^s 

Pm*^m^Ps^s 


(1) 


p ,D  = silicon  and  metal  mass  density 
s m 

respectively 

t_ft  = metal  layer  thickness  and  affected 
m s •* 

semiconductor  depth  respectively 
(cf  Pig.  2) 


Rearranging  Eq.  1 the  ratio  of  the  depth  of  affected  semiconductor 
to  metal  layer  thickness  is: 


P 


<2) 


s 


Hence  the  total  affected  depth  due  to  the  alloying  process  is: 


W 


1 + 


M 


1 - W, 


M 


(3) 


The  eutectic  thickness  at  the  eutectic  temperature  i'  determined 
in  the  same  manner  in  accordance  with  the  weight  percent  at  the 
eutectic  Wg  as: 


(4) 


Therefore  the  regrown  semiconductor  is  just  the  difference  between 
Eqs.  3 and  4: 


(5) 


Since  is  greater  than  Wg  this  is  always  greater  than  zero. 

These  equations  constitute  the  alloy  model  and  allow  calculation 
of  the  various  affected  depths  under  the  given  alloying  conditions 
Another  parameter  of  interest  is  the  density  of  the  eutectic  alloy 
Assuming  a uniform  and  homogeneous  eutectic  region  conservation 
of  mass  gives  the  alloy  density  as: 
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Figure  3.  Solid  Solubilities  of  Impurities 
in  Si  (adapted  from  Ghandi^) 


USE  OF  THE  ALLOY  MODEL 


The  alloy  model  has  been  useful  in  considering  the  types  of 
metal  semiconductor  systems  required  for  an  alloy  junction  target 
Based  on  use  of  the  model,  earlier  writing  experiments  are  better 


understood  enabling  us  to  revamp  our  thinking  on  certain  target 
materials  and  allowing  new  experiments  to  be  outlined. 

One  of  the  first  concepts  for  the  alloy  junction  target  [2] 
was  to  use  a gold  or  silver  metal  layer  on  a silicon  substrate. 

It  was  recognized  that  these  materials  produced  deep  energy  levels 
in  the  silicon  band  gap  and  therefore  would  not  produce  a conven- 
tional diode  as  with  shallow  impurities.  Neverthelessr  it  was  felt 
that  because  of  the  large  weight  percent  of  silicon  in  the  eutectic 
sufficient  penetration  of  the  eutectic  into  the  semiconductor  would 
occur  to  produce  readout  modulation.  Initial  laser  writing  experi- 
ments (3l  verified  that  writing  did  occur  for  a gold  silicon  system, 
however,  readout  modulation  was  poor  when  read  with  the  metal  layer 
on.  At  the  time  of  the  writing  experiments  it  seemed  that  the 

O 

poor  modulation  was  a result  of  too  thick  a gold  film  (1000  A) 
being  used  in  the  experiment  as  well  as  the  grain  structure  thus 
masking  the  effect  of  reading  the  alloyed  region. 

Through  the  use  of  the  alloy  model  it  is  now  believed  that 
there  is  a fundamental  problem  with  systems  such  as  Au-Si  and  Ag- 
Si  that  preclude  their  use  as  alloy  junction  targets.  Based  on 
the  solid  solubility  data  of  Figure  3,  it  is  seen  that  gold  has 
an  extremely  small  solid  solubility  in  silicon  over  the  whole  tem- 
perature range.  What  this  means  in  terms  of  the  target  is  that, 
regardless  of  how  high  above  the  eutectic  temperature  the  system 
is  raised,  there  will  not  be  enough  gold  introduced  in  the  regrowth 
region  to  affect  the  hole  capture  rate.  Therefore,  the  resulting 
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modulation  for  a target  such  as  this  can  only  come  from  the  higher 
electron  stopping  of  the  eutectic  region  relative  to  the  metal 
over  layer.  By  consideration  of  the  model  equations  this  can  be 
shown  to  produce  relatively  weak  modulation.  For  a constant  beam 
energy  (as  in  the  case  of  the  read  beam)  the  mass  range  of  electron 
penetration  is  constant  [4l  and,  therefore; 


PR  = const  = C 

where  p is  the  mass  density  and  R is  the  electron  range. 

Hence  for  the  read  beam  scanning  just  the  metal  layer,  the  range 
to  metal  layer  ratio  is  given  by: 

— = — ^ (8) 

t p t ' 

m ^m  m 

Now  from  Eqs.  6 and  4 it  is  seen  that  the  factor  that  increases 
the  eutectic  thickness  over  the  metal  layer  thickness  shows  up 
as  a reduction  in  the  alloy  density  and  therefore  the  range  to 
eutectic  layer  thickness  is: 


(1  - Wj.) 


or  in  view  of  Eq.  8 


(1  - Wj.) 


(9) 


(10) 
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Therefore  the  reduction  in  modulation  due  to  the  eutectic  is  only 
(1  - Wg)  which  is  insignificant  unless  Wg  becomes  greater  than 
20%  or  so. 

In  view  of  this  result  it  is  now  believed  that  Au  (Wg  = 6%) 
and  Ag  (Wg  = 3%)  would  only  be  useful  if  they  were  loaded  with 
a dopant  material  such  as  Ga  or  Sb.  This  would  allow  the  regrowth 
region  to  be  doped  to  the  opposite  polarity  of  the  top  layer  of 
the  planar  diode  substrate  thus  forming  a surface  diode  with  its 
depletion  field  to  enhance  the  dead  layer  (or  hole  capture)  in 
the  written  area.  The  reason  gold  might  be  used  as  a carrier  for 
the  Ga  or  Sb  rather  than  just  using  the  dopants  themselves  can 
also  be  understood  from  the  model  and  the  phase  diagrams.  Fig- 
ure 4 shows  the  phase  diagrams  for  Ga  - Si  and  Au  - Si.  The  Ga 

o 

- Si  system  not  only  has  a very  low  eutectic  temperature  (=30  C)  , 
but  also  the  slope  of  the  phase  diagram  is  very  steep.  This  means 
very  high  temperatures  must  be  reached  in  order  to  achieve  a deep 
regrown  region.  In  contrast  the  Au  - Si  system  has  a moderately 
high  eutectic  temperature  (370  °C)  with  a much  lower  slope  of  the 
phase  diagram  above  the  eutectic.  The  net  result  is  that  much 
deeper  doped  regions  can  be  achieved  at  lower  temperatures  using 
gold  as  a carrier  for  the  low  solubility  dopants  than  by  using 
the  dopants  themselves.  This  is  illustrated  in  Figure  5 which 
shows  the  regrowth  to  metal  layer  ratiOr  Eq.  5,  as  a function 

of  temperature. 

Experimental  data  are  available  that  not  only  verifies  the 
use  of  gold  as  a dopant  carrier  but  also  verifies  the  alloy  model. 
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Figure  4.  The  Gold-Silicon  and  Gallium  Silicon  Systems 
(From  M.  Hansen  and  A.  Anderko,  Constitution 
of  Binary  Alloys,  McGraw  Hill,  New  York,  1958) 
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Figure  5.  Crystal  Regrowth  for  Alloying  Above 
the  Eutectic  for  Au-Si  and  Ga-Si 

I 

i Figure  6 shows  spreading  resistance  measurements  as  a function 

i 

; of  depth  obtained  from  a diode  formed  under  steady  state  conditions 


I 

* •»- 


using  gold  as  a carrier  for  antimony  dopant.  The  conditions  were 
3.5  mils  of  gold  with  1%  antimony  held  at  a temperature  of  750 
for  one  hour.  Based  on  this  resistance  curve  a eutectic  depth 
of  126  y m and  a regrowth  of  30  ym  is  indicated  giving  a total  af- 
fected depth  of  156  ym.  Calculations  based  on  the  Au  — Si  phase 
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Figure  6.  Spreading  Resistance  on  Depth  Au,  Si  Doped 
- Si  Diode 

diagram  using  the  above  conditions  and  the  model  Equations  (3), 

(4),  and  (5)  give  a eutectic  depth  of  136.02  pm,  regrowth  of  30.23  pm 
and  a total  affected  depth  of  166.25  urn. 

The  model  has  been  used  as  an  aid  in  setting  up  experiments 
to  determine  the  nature  of  the  writing  achieved  under  short  pulse 
conditions.  Based  on  the  excellent  correlation  of  the  model  calcu- 
lations and  the  Au-Sb-Si  diode,  we  know  that  good  abrupt  junction 
diodes  can  be  fojrmed  under  steady  state  conditions.  We  also  know 
from  previous  laser  writing  experiments  that  good,  strong  writing 
effects  can  be  obtained  with  Al-Si  at  temperatures  below  the  melt- 
ing point  of  aluminum  and  with  single  laser  pulses  of  150  urn  sec 
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\ pulse  duration.  Furthermore  from  thermal  model  considerations 

i 

I we  know  that  steady  state  temperatures  can  be  achieved  in  time 

of  the  order  on  nanoseconds.  The  remaining  question  is  whether 

I 

? or  not  the  alloy  process  can  form  good  junctions  within  the  short 

' pulse  times.  Using  the  alloy  model,  we  have  calculated  the  ex- 

pected hole  capture  depth  based  on  two  limiting  cases:  (1)  a diode 

! 

has  been  formed  and  the  capture  depth  is  the  attendant  depletion 
width,  and  (2)  no  diode  is  formed  and  the  captor  depth  is  simply 
the  depth  at  the  regrown  crystal  region.  In  the  case  of  the  latter, 
it  is  assumed  that  the  carrier  lifetime  has  been  killed  in  the 
regrown  region. 

Alloy  junction  diodes  are  usually  modeled  as  an  abrupt  junc- 

I tion  [5]  with  a carrier  concentration  and  space  charge  region  as 

shown  in  Figure  7.  For  this  model  the  depletion  fields  in  the 
n and  p materials  with  no  external  bias  are  given  by  the  equations: 
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q > electronic  charge 

k s Boltzmann  constant 

T « absolute  temperature 

Using  this  diode  model,  the  aluminum  silicon  phase  diagram  of 
Figure  8 and  the  alloy  model,  the  resulting  dead  layer  thickness 
was  calculated  as  a function  of  temperature  for  various  metal 
film  thickness  as  shown  in  Figure  9.  For  these  calculations, 
n on  p planar  diodes  with  constant  n layer  doping  of  10^^/cm^ 
were  assumed.  For  temperatures  just  above  the  eutectic  the  regrown 
layer  of  the  alloy  process  which  provides  the  p doped  portion 
of  the  surface  is  less  than  the  depletion  field  requirement  of 
Eq.  12.  Under  these  conditions  the  p depletion  region  is  just 
equal  to  the  regrown  crystal  region  doped  at  concentration  de- 
termined from  the  solid  solubility  data  of  Figure  3 at  the  spec- 
ified temperature  resulting  in  an  n layer  penetration  of  the  de- 
pletion field  given  by 
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However,  as  the  alloy  temperature  is  raised,  the  regrown  region 
becomes  thicker  than  Xp  of  Eq.  12  and  the  dead  layer  nearly  satu- 
rates with  a slight  slope  from  that  point  on.  This  occurs  in  the 
range  of  585^  to  600  for  the  metal  thicknesses  assumed.  This 
saturation  occurs  because  once  the  regrown  crystal  region  exceeds 
Xp  and  the  full  depletion  width  is  formed,  further  increases  in 
the  dead  layer  come  from  the  increase  in  penetration  of  the  eutectic 
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Figure  8.  The  Aluminum-Silicon  System.  (From  M.  Hansen, 
and  A.  Anderko,  Constitution  of  Binary  Alloys, 
McGraw-Hill,  New  York,  1956.) 

and  regrown  regions  with  increasing  temperatures.  This  is  a slowly 
increasing  function  of  temperature  because  the  regrown  depth  is 
smaller  than  the  depletion  depth  and  increases  in  the  regrown  re- 
gion come  only  from  the  slow  increase  in  the  silicon  solubility 
in  the  Al-Si  eutectic  with  temperature. 
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As  another  possible  bound  for  the  dyanamic  alloy  process  it 
was  assumed  for  calculation  purposes  that  the  dead  layer  was  con- 
fined to  just  the  recrystalized  region.  In  essence  this  assump- 
tion means  that  the  aluminum  damages  the  regrowth  region  but  does 
not  go  in  substitutionally  to  form  an  active  dopant  ion.  This 
results  in  a dead  layer  of  tg  from  Eg.  2.  Results  of  calculations 
based  on  this  model  are  shown  in  Figure  10. 

Based  on  these  assumptions,  there  are  significant  differences 
between  calculated  dead  layers  for  metal  layer  thicknesses  of  250, 

O 

500,  and  1000  A and  these  differences  should  be  experimentally 
visible.  In  view  of  this,  laser  writing  experiments  with  a newly 
developed  single  pulse  mode  for  more  accurate  temperature  rise 
calibration  have  been  initiated.  These  experiments  and  their  re- 
sults will  be  discussed  in  the  next  reporting  period.  By  writing 
at  various  power  settings  (i.e.  varying  temperature)  and  with  dif- 
ferent film  thickness,  subsequent  electron  probe  analysis  of  the 
dead  layer  should  give  the  nature  of  the  short  pulse  alloy  process 
when  compared  with  the  alloy  model. 

SUMMARY 


An  alloy  model  has  been  developed  which  describes  the  steady 
state  alloying  process  allowing  calculation  of  the  various  crystal 
regrowth  and  eutectic  penetrations.  Application  of  the  model  to 
an  Au-Sb-Si  diode  verifies  its  correctness.  The  model  has  been 
used  to  aid  in  the  understanding  of  basic  metal  semiconductor  re- 
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quirements  for  an  alloy  junction  target.  From  the  alloy  point 
of  view  these  requirements  translate  into  a low  temperature  eutectic 
in  conjunction  with: 
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1)  A high  weight  percent  of  Si  in  the  eutectic  (with  a 
minimum  of  at  least  5%) 

2)  High  solid  solubility  of  the  metal  or  a carried  dopant 
in  order  to  (a)  dope  the  crystal  (b)  cause  sufficient 
damage  in  the  regrowth  region  to  produce  a deep  dead  layer. 


^ SECTION  III 
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\ THERMAL  EQUILIBRIUM  STUDIES 
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During  this  reporting  period,  thermal  equilibrium  studies 
of  the  alloy  junction  process  were  carried  out.  Planar  diode 
structures  were  metalized  in  30  mil  dot  patterns  and  the  devices 
were  examined  in  a scanning  electron  microscope  using  both  second- 
ary and  diode  signal  readout,  before  and  after  heating  in  a tem- 
perature regulated  furnace.  The  dot  patterns  were  formed  by  evap- 
oration or  sputtering  of  thin  metal  layers  through  a metal  mask. 
Both  high  and  low  temperature  regimes  were  investigated. 

Since  the  gold  eutectic  with  silicon  is  a system  of  consider- 
able interest,  thin  films  of  this  metal  on  archival  diode  sub- 
strates were  studied.  Figure  11  illustrates  a typical  experiment 
consisting  of  several  dot  arrays  produced  by  sputter  depositing 
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METALIZED  DOTS. 
Si  UNDERNEATH 
SPUHER  ETCHED 

METALIZED  DOTS 


04^ — CONTACT 


BEFORE  HEATING 


AFTER  HEATING 


Figure  11.  Illustration  of  Typical  Thermal  Equilibrium 
Alloy  Junction  Experiment,  Consisting  of 
Small  Metalized  Dot  Arrays  Deposited  on 
Sputter  Etched  and  Unetched  Regions  of 
a Planar  Diode  Device 

200  A of  gold  on  PAR  lA  target.  PAR  lA  targets  were  produced  by 
two  separate  boron  implants  into  an  n-type  substrate  to  form  a 
p^pn  structure.  In  these  experiments,  all  substrate  material  was 
<111>,  since  the  solid-liquid  interface  lies  in  this  plane. 
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A high  degree  of  surface  uniformity  is  required  for  thin, 
uniform  regrowth  in  the  alloyed  region  Even  etch  cleaning  of 
already  polished  wafers  can  cause  sufficient  surface  irregularity 
to  result  in  nonuniform  junction  depths.  To  avoid  this  spiking 
action,  fine  polishing  (0.1  y)  compounds  were  used  in  the  production 
of  alloy  junction  devices  to  prepare  the  surface  for  metalization. 
This  step,  however,  can  result  in  damage  depths  =ly,  which  would 
create  a "dead  layer"  perceptible  during  readout  of  the  device 
with  an  electron  beam. 

Sputter  cleaning,  therefore,  was  used  to  remove  any  native 
oxide  on  the  structure  surface.  This  method  was  employed  to  im- 
prove wetting  of  the  metal  and  enhance  alloy  formation  as  well. 

Both  sputter  etched  and  unetched  areas  were  metalized  for  compar- 
ison. Diode  readout  measurements  were  taken  for  both  types  of 
areas  and  the  bar  diode.  No  distinguishable  dead  layer  was  observed 
for  the  etched  versus  unetched  areas,  as  shown  in  Figures  12  and 
13. 


Sample  heating  was  performed  in  a furnace  in  a quartz  boat 
located  in  the  flat  zone,  surrounded  by  argon  ambient.  Target 
structures  were  maintained  at  temperature  5 minutes  before  removal. 
For  sample  PAR  lA,  the  low  temperature  run  was  performed  at  375  °C, 
5 °C  above  the  eutectic  temperature.  Figures  12  and  13  compare 
the  diode  signal  before  and  after  baking  for  the  metalized  dot 
arrays,  and  Figure  14  indicates  the  signal  before  and  after  heat- 
ing in  the  untreated  diode  region.  Little  evidence  for  alloying 


was  observed. 
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READING  ELECTRON  BEAM  ENERGY  (k«V) 

Figure  12.  Collection  Efficiency  from  P^pn  Planar  Diode 
Device  Before  and  After  375  °C  Baker  in  a 
Region  Metalized  with  200  A of  Gold,  Sputter 
Cleaned  Before  Deposition 
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Figure  13.  Collection  Efficiency  from  p^pn  Diode  Device 
Before  and  After  375  *^C  Bake#  in  a Region 
Metalized  with  200  A of  Gold. 
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Figure  14.  Collection  Efficiency  from  a p^pn  Diode 
Before  and  After  375  °C  Bake 

When  this  experiment  was  repeated  for  the  375  °C  bake  tempera- 
ture, agglomeration  of  the  gold  was  observed,  particularly  in  the 
sputter  etched  region.  Figures  15,  16,  17,  and  18  illustrate 
secondary  electron  micrographs  taken  before  and  after  bake  showing 
changes  in  the  gold  metalization  before  and  after  375  °C  bake, 
in  sputter  etched  and  unetched  regions.  Before  the  bake  cycle, 
the  200  A sputter  deposited  gold  film  is  smooth.  The  rim  of  the 
sputter  etched  region  is  the  only  visible  structure.  After  bake, 
the  gold  film  is  no  longer  uniform,  particularly  in  the  dot  de- 
posited on  a sputter  etched  region.  In  the  non-sputter  etched 
region,  the  gold  appeared  to  pull  away  from  regions  with  poor 
wetting  properties.  Figure  9 shows  a close-up  of  the  sputter 
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Figure  15.  Secondary  Electron  Micrograph  Showing  Dot  A4 
on  Target  PAR  1(11) -2  Before  375  °C  Bake. 
Surface  was  sputter  etched  before  deposition 
of  200  A gold  film. 


Figure  16.  Secondary  Electron  Micrograph  Showing  Dot  B4 
on  Target  PAR  1(11) -2  Before  375  °C  Bake. 
Surface  was  not  sputter  etched  before 
deposition  of  200  A gold  film. 
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Figure  17.  Secondary  Electron  Micrograph  Showing  Dot  A4 
on  Target  PAR  1(11) -2  after  375  ”c  bake. 


Figure  18.  Secondary  Electron  Micrograph  Showing  Dot 
on  Target  PAR  1(11) -2  after  375  °C  bake 
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Figure  19.  Secondary  Electron  Micrograph  of  Faceted  Gold 
Islands  Present  in  Sputter  Etched  Region  of 
Target  PAR  1(11) -2  After  375  Bake.  The 
gold  regions  are  several  microns  in  diameter. 

etched  region  which  displays  clumps  of  gold  across  the  surface. 
Multifaceted  crystalline  gold  areas  are  present  over  the  entire 
region.  These  islands  are  several  microns  in  diameter.  This  ef- 
fect could  be  the  combined  result  of  surface  tension  forces  and 
the  realignment  of  the  gold  film  with  the  intermediate  silicide 
layer  whose  growth  is  accelerated  by  the  heating  process.  This 
growth  rate  is  also  dependent  on  surface  conditions.  An  unmetal- 
ized  device  which  had  been  sputter  cleaned  was  carefully  examined 
using  the  secondary  signal  in  a scanning  electron  microscope  to 
determine  if  the  process  induced  surface  irregularities.  Even 
at  the  highest  magnifications,  no  such  nonuniformities  were  ob- 
served. 
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High  temperature  (800  °C)  bakes  were  also  performed  to  attempt 
verification  of  the  thermal  equilibrium  model  (described  in  Sec- 
tion II) . Again,  secondary  and  diode  measurements  were  carried 
out  before  and  after  heating  of  1000  A gold  films  on  target  Par  1 
(11) -4.  Again,  after  the  bake,  the  gold  observed  to  be  smooth 

before  heat  treatment,  agglomerated.  Regular  triangular  gold  crystals 
were  visible  in  the  region  sputter  etched  before  gold  deposition, 
whereas  in  the  nonsputter  etched  region,  the  gold  had  simply  coagu- 
lated into  blobs  (Figures  20  - 25) . Diode  measurements  were  made 
on  the  gold  islands;  however,  uncertainty  in  the  thickness  of  these 
faceted  islands  made  a direct  determination  of  the  affected  depth 
difficult.  There  was  evidence  that  alloying  did  occur  underneath 
the  coagulated  gold  regions;  however,  the  data  was  of  insufficient 
quality  to  allow  corroboration  with  the  model. 


Figure  20.  Secondary  Electron  Micrograph  of  Dot  B6  of 

Target  PAR  1(11) -4,  Formed  by  Sputter  Deposition 
of  1000  A of  Gold,  Before  800  Bake 
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Figure  21.  Secondary  Electron  Micrograph  of  Dot  AS  of 

Target  PAR  1(11) -4,  Formed  by  Sputter  Deposition 
of  1000  A of  Gold  in  a Sputter  Etched  Region, 
Before  800  °C  Bake 


Figure  22.  Secondary  Electron  Micrograph  of  Dot  B6  of 
Target  PAR  1(11) -4  After  800  °C  Bake  for 
5 Minutes  in  Argon  Ambient,  Showing 
Agglomeration  of  Gold 


Figure  23.  Secondary  Electron  Micrograph  Giving  Magnified 
Views  of  Regions  of  Agglomerated  Gold  in  a 
Non-sputter  Etched  Region  of  Target  PAR  1(11) -4, 
Dot  B6 


Figure  24.  Secondary  Electron  Micrograph  of  Dot  B6 
on  Target  PAR  1(11) -4  After  800  C Bake 
’ for  5 Minutes  in  Argon,  Showing  Triangular 

Platelets  of  Gold 
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Figure  25.  Secondary  Electron  Micrographs  Giving 
Magnified  Views  of  Triangular  Gold 
Platelets  Formed  in  Sputter  Etched 
Regions  of  Target  PAR  1(11) -4,  Dot  A5. 
Platelets  are  typically  2 - 3 p in 
diameter . 

Aluminum  dot  arrays  on  silicon  planar  diode  structures  were 
also  investigated,  since  this  eutectic  system  is  also  of  consider- 
able interest.  For  these  experiments,  aluminum  was  sputtered  in 
a dot  pattern  on  n and  p diode  structures.  When  this  1000  A alu- 
minum film  was  examined  at  high  magnification  and  low  beam  volt- 
ages in  the  scanning  electron  microscope,  collection  efficiencies 
resulted  which  indicated  the  presence  of  tiny  pinholes.  While 
too  small  to  be  observed  with  secondary  electrons,  they  did  allow 
penetration  of  the  beam.  The  presence  of  such  pinholes  may  be 
desirable,  since  they  may  reduce  the  thermal  conductivity  in  the 
film.  Similar  effects  had  been  seen  in  gold  sputter  deposited 
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Evaporated  films  were  also  used  in  these  studies.  A 1000  A 
thick  aluminum  film  was  deposited  on  target  MP  3-4(3)  for  investi- 
gation of  the  affected  alloy  depth  at  790  °C  heating  in  argon  am- 
bient for  5 minutes.  Again,  diode  and  secondary  measurements  were 
made  before  and  after  the  furnace  cycle. 

These  measurements  indicated  that  a diode  had  indeed  been 
produced  by  alloy  formation,  but  again,  agglomeration  of  the  film 
during  bake  resulted  in  difficulty  in  obtaining  absolute  measure- 
ments of  the  alloy  depth.  These  agglomerations  were  in  the  form 
of  film  discontinuity  and  triangular  platelets  of  the  order  of 
10 y m in  size.  No  sputter  cleaning  had  been  done  on  this  sample. 
Secondary  electron  micrographs  illustrating  these  metal  islands 
are  shown  in  Figures  26  and  27.  Figure  28  shows  diode  signal  from 


Figure  26.  Secondary  Electron  Micrograph 
Showing  Restructured  1000  A 
Aluminum  Film  on  Target  MP  3-4(3) 
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Figure  27.  Secondary  Electron  Micrograph  Showing 
Magnified  A1  Triangular  Platelet 
( ~ 20vj  on  a side)  from  Target 
MP  3-4(3)  After  Heating  at  690  °C 


Figure  28.  Diode  Readout  from  a Region  Heated 
at  690  °C  on  Target  MP  3-4(3), 
Showing  Low  Gain  in  Regions  of  Metal 
Agglomeration 
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a similar  region.  The  darker  the  signal,  the  lower  the  gain  in 
an  area,  such  as  at  the  location  of  the  triangles. 

SUMMARY 

The  experiments  did  clearly  show  that  writing  is  not  accom- 
plished by  simply  bringing  the  constituents  to  their  eutectic 
temperature.  Higher  temperatures  are  required;  but  the  experi- 
ments neither  proved  nor  disproved  the  predictions  of  the  model 
in  Section  II  on  this  point. 

Experiments  did  indicate  that  thermal  biasing  of  this  type 
of  target  would  be  quite  difficult,  in  view  of  the  formation  of 
silicide  layers,  and  of  the  restructuring  of  the  overlaying  films. 
In  addition,  the  experiments  showed  that  there  is  a definite  prob- 
lem with  making  the  metal  species  wet  the  surface  of  the  silicon 
during  alloying.  For  the  alloy  function  approach  to  work,  a method 
must  be  found  to  make  the  metal  stay  put  when  heated  so  that  alloy- 
ing can  take  place.  The  interaction  of  the  metal  layers  with  de- 
fects during  the  alloy  process  is  another  potential  problem. 


SECTION  IV 


THERMAL  MILLING  OF  Se-As  ALLOYS 


INTRODUCTION 


The  electron  milling  concept  is  shown  schematically  in  Fig- 
ure 29.  A focused  electron  beam  locally  heats  a film  of  Se  on 
a Si  diode  target  causing  vaporization  of  the  Se.  After  writing, 
the  target  is  read  using  an  electron  beam  of  the  same  energy 
but  of  lower  beam  current.  Regions  of  reduced  Se  thickness  due 
to  the  machining  produce  higher  readout  signals  from  the  diode 
because  of  the  increased  penetration  of  the  electron  beam  energy 
in  the  Si  diode. 


WRITE  READ 

Figure  29.  Electron  Beam  Milling  Concept 


Se  alloyed  with  As  was  chosen  for  this  application  because 
it  is  a material  with  the  following  properties:  low  thermal 
conductivity,  moderate  mass  density,  low  boiling  temperature, 
good  stability,  and  uniform  amorphous  thin  films  are  possible. 
Amorphous  Se  films  are  known  to  have  a strong  tendency  to  recrys- 
tallize at  room  temperature.  Crystalline  films  are  in  general 
undesirable  because  they  have  higher  thermal  conductivities  and 
tend  to  be  nonuniform  on  the  scale  of  0.1  urn.  Also  films  which 
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undergo  amorphous  to  crystalline  transitions  near  room  temperature 
are  unlikely  to  be  dimensionally  stable  over  long  periods  of 
time.  The  addition  of  As  to  Se  has  been  known  for  some  time 
to  significantly  stabilize  Se  against  crystallization. 

CHOICE  OF  SUBSTRATE  PLANAR  DIODE 


N type  planar  diodes  were  chosen  as  the  substrates  for  the 
Se  milling  investigations.  These  are  identical  to  the  As^  implant 
through  oxide  targets  developed  during  Phase  I and  Phase  II.  They 
are  processed  using  the  NAR  process  described  in  the  Phase  II 
Interim  Report  (June  1 - December  31,  1977) . These  diodes  have 
high  electron  beam  gain  at  low  beam  energies  and  are  suitable  for 
general  writing  studies  as  well  as  substrates  for  high  resolution 
recording . 

PREPARATION  OF  Se-As  ALLOY  FILMS 


The  optimum  choice  of  the  Se-As  alloy  composition  for  this 
application  has  not  been  determined.  In  general,  the  thermal 
conductivity  of  Se-As  glasses  increases  with  progressively  greater 
As  content  up  to  40%  As  ( 6 J.  This  effect  is  less  than  a factor 
of  two  so  that  other  factors  will  dominate  the  choice.  In  particu- 
lar, the  conditions  of  film  preparation  may  be  more  important. 

As  mentioned  above,  the  As  stabilizes  the  amorphous  Se  films 
against  recrystallization  (14,  17,  18,  19,  20,  33].  Amorphous 
Se  exists  in  both  Seg  rings  and  long  chains.  Both  of  these  forms 
tend  to  form  complex  crystalline  structures  with  higher  thermal 
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conductivities.  Monoclinic  and  hexagonal  crystals  are  formed 
by  the  Se  rings  and  long  chains,  respectively.  Addition  of  As 
forms  AsSe2^2  that  break  up  the  rings  and  cause  cross- 

linking  of  the  chain.  1131.  There  is  theoretical  (201  and  ex- 
perimental evidence  [331  that  only  100  ppm  of  As  is  required  to 
achieve  this  stabilization.  At  higher  As  concentrations  this 
process  eventually  results  in  the  As2Se2  structure.  Heating 
above  the  glass  transition  temperature  Tg  may  cause  some  crys- 
tallization (6,  24  1 in  the  higher  As  concentration  alloys.  The 
extent  to  which  this  crystallization  persists  when  subjected  to 
rapid  or  slow  cooling  is  not  known. 

Based  on  this  evaluation  it  was  decided  to  produce  As-Se 
alloys  with  compositions  ~20%  As.  The  principal  reasons  for 
this  were:  this  is  the  eutectic  composition  and  so  it  has  the 
lowest  melting  point,  and  it  was  the  composition  utilized  in 
the  previous  GE  work  at  Electronics  Laboratory,  Syracuse,  N.Y. 

As  will  be  seen  below,  all  of  the  films  prepared  to  date  have 
much  lower  As  concentrations  (~1%).  These  films  seem  to  be 
completely  suitable  as  memory  materials.  As  just  discussed, 
there  is  evidence  that  As  concentrations  as  low  as  100  ppm  may 
produce  equally  or  more  stable  amorphous  Se  films  [19,  20 ] . 

One  percent  As  films  are  apparently  a feasible  composition  for 
xerography  drums  (12,  111  . with  a heating  rate  of  3 °C/min  the 
crystallization  temperature  of  As-Se  increases  from  160  to 
196  °C  as  the  As  concentration  is  increased  from  0 to  2%  (20  ]. 
This  process  provides  some  indication  of  the  strong  effect  of 
As  on  the  crystallization. 


As-Se  alloys  were  prepared  by  heating  a charge  of  23.3  g 
of  Se  and  5.8  g of  As  in  an  open  carbon  crucible  in  flowing  argon 
at  atmospheric  pressure.  For  the  first  alloy  run  the  temperature 
was  increased  slowly  over  1 hour  to  315  °C  and  then  dropped  back 
to  300  °C  for  12  hours.  Considerable  evaporation  of  material 
occurred.  Subsequent  X-ray  florescence  analysis  of  the  alloy 
indicated  a composition  of  1%  As  with  no  other  impurities  above 
0.1%.  Electron  microprobe  studies  indicated  no  composition  varia- 
tions in  the  alloy  although  the  1%  As  composition  is  below  the 
detection  threshold  due  to  the  Bremsthrahlung  background.  It 
is  believed  that  the  As  preferentially  evaporated  from  the  cru- 
cible. The  vapor  pressure  of  As  in  a dilute  Se  alloy  is  not 
known  but  it  is  predicted  that  the  As  and  Se  vapor  pressures 
should  be  equal  at  387  °C  [12].  Hence  it  is  possible,  but  un- 
likely, that  the  As  would  preferentially  evaporate  at  300  °C. 

A second  possibility  is  the  formation  of  As  oxides  due  to  an 
impure  Ar  gas  or  other  leakage  in  the  system.  All  of  the  As 
oxides  are  unstable  or  sublime  below  315  °C.  A second  alloy 
run  at  maximum  temperature  of  230  °C  did  not  produce  mixing  of 
the  As  and  Se.  Future  alloys  will  be  made  in  sealed  evacuated 
quartz  ampoules  at  temperatures  - 350  °C  to  prevent  loss  of  As 
during  alloying.  All  of  the  Se  films  prepared  for  the  investiga- 
tion described  below  used  this  material  from  the  first  alloy 
run. 


The  films  were  evaporated  in  a bell  jar  using  indirectly 
heated  carbon  crucibles.  The  thickness  was  monitored  with  an 
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Infocon  quartz  crystal  monitor  and  the  evaporation  controlled 

O 

by  a shutter.  Deposition  rates  were  typically  10  to  100  A/sec 
at  the  target.  The  films  have  not  yet  been  analyzed,  but  based 
on  the  work  of  Sigai  [12]  and  our  estimates  of  the  boat  tempera- 
ture {-  300  °C) , they  should  be  ^ 1%  As.  The  Se  should  preferen- 
tially evaporate  because  of  its  higher  vapor  pressure.  Since 
the  shutter  is  in  place  during  the  initial  part  of  the  evaporation 
some  As  enrichment  of  the  charge  occurs  before  evaporation  is 
started.  Sigai  found  that  for  dilute  (0.5%  As)  alloys  it  was 
impossible  to  reach  a temperature  high  enough  to  obtain  equal 
evaporation  rates  for  As  and  Se.  This  he  attributed  to  the  poor 
thermal  conductivity  of  the  Se  and  the  effect  of  evaporative 
cooling.  Because  of  the  compensating  enrichment  effect  and  the 
fact  that  the  charge  was  typically  evaporated  to  near  completion 
it  is  expected  that  the  films  have  ~ 1%  As.  Nicroprobe  dispersive 
analysis  indicated  less  than  2%  As  in  the  two  samples  tested. 
(NARl-2(2),  NAR1-12(1)  ). 


Test  glass  disks  are  included  with  every  evaporation.  These 
are  then  covered  with  a thin  A1  film  and  the  physical  thickness 
measured  using  a Sloan  M-lOO  Angstrometer  (Fizeau  technique) . 
Figure  30  shows  the  correlation  between  the  reported  thickness 
based  on  the  quartz  monitor  and  the  physical  thickness.  The 

2 

quartz  monitor  thickness  is  based  on  a mass  density  of  4.25  g/cm 

2 

[15]  and  includes  cos e and  r corrections  for  the  differences 
in  monitor  and  sample  positions.  The  slight  disagreement  is 
most  likely  due  to  a non  cosine  evaporation  distribution  for 
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the  carbon  crucible  since  about  10%  of  the  Se  charge  is  shielded 
from  the  monitor  by  the  edge  of  the  crucible. 


OPTICAL  THICKNESS  1^) 


Figure  30.  Comparison  of  Evaporated  Se  Film  Thick 
hess  Measured  by  Quartz  Monitor  and 
Optical  Interference 


STABILITY  OF  THE  EVAPORATED  FILMS 
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Optical  and  SEM  examination  of  the  evaporated  films  indicate 
that  they  are  amorphous  and  of  generally  good  quality.  Test 
films  evaporated  on  glass  disks  and  stored  in  a normal  laboratory 
environment  for  several  months  still  exhibit  the  reddish  color 
and  good  optical  transmission  characteristic  of  amorphous  Se 


films.  No  degradation  of  the  films  has  been  observed  although 
those  films  which  had  A1  evaporated  on  them  have  shown  some  re- 
action with  the  Al.  A1  is  a very  reactive  material  and  so  this 
is  not  surprising. 

Se  WRITING  EXPERIMENTS 


Electron  beam  writing  has  been  done  in  the  Coates  and  Welter 
(C&W)  field  emission  electron  writing  column  purchased  during 
Phase  I.  For  writing,  the  beam  is  stepped  in  a 16  x 16  matrix 
of  dots  using  two  synchronized  staircase  generators.  A number 
of  modes  of  operation  are  possible.  The  beam  moves  between  the 
dots  of  the  matrix  in  <1  usee.  The  time  between  transitions 
to  a new  spot  can  be  varied  between  20  msec  and  0.088  msec. 

Without  the  use  of  beam  blanking,  the  dwell  time  per  spot  can 
be  varied  between  these  limits.  The  matrix  generator  can  be 
synchronized  to  a pulse  generator  which  is  used  to  drive  the 
electrostatic  beam  blanker  of  the  column.  In  this  mode  the  beam 
dwell  time  can  be  as  small  as  0.5  gsec.  This  minimum  time  is 
limited  by  the  performance  of  the  C6W  blanker.  The  beam  blanker 
reduces  the  beam  current  to  about  10%  which  is  normally  sufficient 
to  prevent  further  writing.  In  either  mode,  single  lines  or 
the  entire  matrix  can  be  written  in  a single  shot  or  continuous 
modes. 


In  general,  the  operation  of  the  column  has  been  erratic. 
This  has  severely  hampered  progress  on  electron  milling.  The 
typical  best  performance  of  the  column  is  100  nA  at  15  keV  with 


a spot  size  ]f0.1ym  and  25  nA  at  5 keV  with  a similar  spot  size 
> 

~0.1y.  At  higher  currents  the  spot  size  increases  draunatically . 
The  introduction  of  the  beam  blanking  aperture  into  the  column 
tends  to  further  degrade  the  beam  current  and  spot  size. 


l 
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Initially  attempts  were  made  to  write  on  the  0.05  wm  Se 
film.  These  were  generally  unsuccessful  although  some  writing 
was  obtained.  Figure  31  shows  a secondary  image  of  a portion 
of  the  16  X 16  dot  array  written  with  70  nA  at  5 keV  and  a dwell 
time  of  20  msec/spot  repeated  4 times  at  each  spot.  The  spot 
size  was~0.2  to  0.3pduring  this  writing  experiment.  Each  of 
the  3 white  dots  corresponds  to  a single  written  bit.  Estimates 
of  the  temperature  rise  achieved  in  this  case  are  200  °C  which 
is  just  sufficient  to  melt  the  Se  film.  The  electron  beam  po- 
sition is  not  completely  stable  because  of  a small  differenti- 
ated part  of  the  TV  raster  sync  pulse  signal  which  couples  into 
the  deflection  coils  of  the  column.  This  couples  equally  into 
each  axis,  so  for  long  dwell  times  the  beam  spends  3%  of  the 
time  deflected  at  a 45°  angle.  This  causes  a slight  heating  of 
the  film  in  these  adjacent  areas.  Apparently  because  of  the 
surface  tension  character itics  of  the  molten  Se,  it  is  pulled 
towards  these  slightly  warmer  areas  resulting  in  two  raised  regions 
of  Se  on  each  side  of  the  centrally  heated  spot.  The  origin 
of  the  third  raised  spot  is  not  as  clear  but  may  be  related  to 
astigmatism  in  the  spot.  Using  diode  signals  it  is  inferred 
that  the  amount  of  change  in  Se  thickness  is  ^ 5%. 
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Figure  31.  Secondary  Emission  Image 
of  Viscous  Flow  Writing 
on  a 0.05  Micron  Thick 
As-Se  Film.  Electron 
beam  writing  at  5 keV  and 
70  nA  with  a spot  size 
~0.25  microns. 


This  does  demonstrate  ~0.3  ^writing  of  Se  films.  The  time 
required  to  form  these  structures  is  20  msec.  At  shorter  times 
no  writing  was  seen.  Several  pulses  at  20  msec  were  required 
to  produce  a visible  effect.  Se  is  a very  viscous  material  ( 16  ] 
and  hence  flows  very  slowly  at  these  temperatures.  This  type 
of  viscous  flow  writing  is  not  useful  for  memory  applications. 


In  its  present  condition  the  column  does  not  provide  suf- 
ficient power  to  vaporize  thin  Se  films.  In  order  to  study  vapor- 
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Figure  32.  (a)  Secondary  Emis- 

sion Image  of  a Dot 
Matrix  Written  on 
1.8  Micron  Thick 
As-Se  Film  at  100  nA, 
16  keV  and  100  ysec. 
per  spot. 


(b)  Typical  Diode  Signal  Readout 
of  a Written  Dot  Matrix.  Right 
trace  is  beam  on  signal  in  an 
area  with  no  Se  film,  middle 
trace  is  the  modulated  diode 
signal  as  the  beam  is  scanned 
across  a row  of  dots.  The  base 
line  trace  is  the  beam  off  refer 
ence.  Left  image  is  a diode  sig 
nal,  Z-axis  modulated  image,  of 
the  written  dot  matrix. 


ization  writing  it  was  necessary  to  study  thicker  Se  films  and 
to  raise  the  beam  voltage.  Much  greater  temperature  rises  cat: 
be  achieved  in  thicker  films  because  of  the  increased  distance 
to  the  high  thermal  conductivity  silicon  substrate.  Figure  32 
shows  a secondary  emission  image  of  some  typical  writing  on  the 
1.8  Mm  Se  films  written  at  15  keV  and  140  nA.  The  electron  spot 
size  was  0.2  pm.  The  limitation  on  the  size  of  the  written  regions 
is  ~ ] M (3ue  to  electron  scattering. 
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The  depth  of  the  written  holes  is  determined  from  the  electron 
beam  induced  current  in  the  underlying  pn  junction.  This  provides 
a very  convenient  method  for  the  determination  of  the  residual 
Se  at  the  bottom  of  the  hole.  Such  measurements  would  be  nearly 
impossible  by  other  techniques.  Figure  33  shows  a typical  set 
of  such  measurements.  Shown  is  the  diode  collection  efficiency 
in  the  center  of  the  written  dots  measured  at  several  beam  voltages. 
For  writing  the  electron  beam  voltage  was  15  keV  and  140  nA  or 
70  nA.  The  spot  size  was  ~0.1m  m so  the  heated  volume  is  de- 
termined by  the  electron  range.  The  diode  measurements  are  done 
at  low  beam  currents  ~10  nA  to  prevent  further  machining  of  the 
film.  No  degradation  or  change  in  signal  were  seen  during  these 
measurements.  The  solid  curves  are  a model  calculation  based 
on  the  following  assumptions.  The  collection  probability  for 
the  diode  is  as  shown  in  Figure  19  of  the  Phase  I Fourth  Interim 
Report  (January  16  - April  15,  1977)  for  Sv*10®  cm/sec.  This 
is  appropriate  for  the  NAR  type  diodes  used  here  and  provides 
a good  fit  to  the  measured  collection  efficiencies  as  shown  in 
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Figure  33.  Diode  Gain  Measurements  for  the  Determination 
of  Residual  Se  Film  Thickness  in  Milled  Holes 
for  Various  Indicated  Beam  Currents  and  Dwell 
Times ■0.24  msec,  and  140  nA,  • 16  msec,  and 
140  nA  (X)  1 msec,  and  140  nA,  (0)  0.11  msec 
and  140  nA,  (A)  16  msec,  and  70  nA  and  (0) 

1 msec,  and  70  nA.  All  writing  at  16  keV. 

Model  curves  based  on  Equation  (14). 

Figure  17  of  the  same  report.  The  effect  of  the  Se  is  to  absorb 
electron  beam  energy  as  though  it  were  a silicon  film  of  thickness 

(pSe'^PSi^ "^Se  PSe  PSi  mass  densities  of  Se  and 

Si  respectively,  and  Tg^  is  the  residual  Se  thickness  after  machin- 


ing. The  electron  energy  loss  is  described  by  the  Grun  formulation 
as  described  on  page  100  of  the  Phase  I First  Quarterly  Report 
(April  15  - July  15,  1976).  This  is  a reasonable  although  not 
exact  method  of  adapting  the  uniform  media  Grun  function  data 
to  this  two  layer  case.  It  is  identical  to  the  assumptions  used 


in  the  thermal  modeling  calculations  described  in  more  detail 
in  Section  V of  this  report.  Using  these  assumptions  and  gener- 
alizing to  two  films  of  thickness  and  T2  with  mass  densities 
and  P2  the  collection  efficiency  is  given  by: 

1.1 

^ ^'‘3)  P{q-qj-q2)RG(psi'EB)  (14) 

qi+q2 


q2  = - -2 


Rg(C'2  ^b) 


where  R^(p,Egj  ^ ^ E^  in  KeV  6 = 4.28xl0~®  g/cm^  (15) 

is  the  Grun  range  which  depends  upon  the  mass  density  and  the 
beam  energy  Eg.  X (q)  is  the  normalized  Grun  function  for  which 
we  use  the  computational  approximation  [36) 

X (q)  = 0.6  + 6.21q  - 12.40q^  + 5.69q^  (16) 

As  can  be  seen  in  Figure  33,  the  fit  to  the  data  is  satisfac- 
tory. The  systematic  deviations  are  piohiably  due  to  the  error 
in  approx  im.at  ing  the  energy  loss  function.  These  problems  notwith- 
standing, the  data  in  Figure  33  still  provide  a useful  measure 
of  the  residual  Se  left  in  the  written  areas.  If  the  thickness 


is  always  cstirT.ated  for  approximately  the  same  collection  efficiency 
(••0.1)  the  thickness  determinations  will  be  the  most  reliable 


and  have  the  least  scatter.  Table  1 shows  the  measured  writing 
depths  for  three  different  beam  currents  as  determined  from  Fig- 
ure 33. 
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TABLE  1 

MEASURED  WRITING  DEPTHS  FROM  FIGURE  5 
FOR  16  keV  WRITING  ON  1.6 u As-Se 


Beam 

Dwell 

Residual 

Current 

Time 

Thickness 

yA 

Visec 

y 

140 

16 

0.12 

140 

1 

0.32 

140 

0.24 

0.34 

140 

0.11 

0.45 

70 

16 

0.28 

70 

1 

0.45 

Secondary  emission 

images 

of  additional  writing  experiments 

show  that  the  writing  depth  becomes  very  small  for  dwell 

times 

between  10  and  50  u sec 

MODELS  FOR  THE  WRITING 

for  beam 

PROCESS 

currents  in  the  range  of 

100  nA. 

' I 
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A number  of  factors  which  are  important  for  an  understanding 
of  the  thermal  milling  process  will  be  considered  separately. 

These  are;  the  vapor  pressure  and  its  relation  to  the  evapora- 
tion rate  of  the  material;  heating  models  for  the  temperature 
rise  in  the  multilayer  targets;  the  effect  of  the  latent  heats 
of  vaporization  and  fusion;  and  the  thermal  properties  of  amorphous 
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Se-As  alloys.  Baseo;  on  these  discussions  a quantitative  model 
for  the  milling  process  will  be  described. 

A.  Thermal  Properties 

The  thermal  conductivity  of  bulk  amorphous  Se-As  alloys 
is  summarized  in  Table  2 based  on  the  measurements  of  several 
investigators.  As  can  be  seen  from  the  table,  there  is  consider- 
able scatter  in  the  values  obtained.  Some  of  the  scatter  is 
undoubtedly  due  to  the  fact  that  Se  can  have  many  amorphous  struc- 
tures. (See  Reference  20.)  Also  note  that  all  of  these  measure- 
ments are  on  bulk  material,  which  is  prepared  by  quenching  or 
slow  cooling  from  the  melt  rather  than  vacuum  evaporation.  For 
our  modeling  studies  a value  of  0.005  watt/cm/sec  will  be  adopted 
as  a nominal  pessimistic  value.  The  thermal  conductivity  of 
Se-As  tends  to  increase  moderately  with  temperature  except  for 
a dip  at  the  glass  temperature  Tg  which  occurs  between  50  and 
200  °C  depending  upon  the  composition  (61 . 

The  heat  capacity  is  of  less  critical  importance  since  it 
determines  only  the  thermal  time  constant  and  not  the  maximum 
temperature  rise.  A good  nominal  value  at  room  temperature  for 
amorphous  Se  is  Cp  * 1.35  J/cm^/deg  C (221. 

B.  Vapor  Pressure 

The  vapor  pressure  curves  (24,  23,  121  for  pure  As  and  Se 
are  shown  in  Figure  34.  The  standard  vapor  pressure  curve  for 
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As  is  based  on  the  sublimination  of  As  from  the  solid  phase. 

Sigai  [12]  has  estimated  from  thermodynamic  data  the  correction 
for  the  vapor  pressure  over  the  liquid  phase.  This  correction  is 
very  small.  The  Se  curve  can  be  described  by 

log  * A - B/T  (17) 

where  denotes  the  vapor  pressure  in  units  of  atmospheres, 

A * 5.6  and  B = 5376  °K. 


\i 


For  significant  concentrations  of  As,  the  vapor  pressure 
of  each  component  of  the  alloy  can  be  determined  ng  the  law 
of  Raoult: 


P - P 
*^As  *^As 


= 1-X 


P - P 
Se  *^Se 


As 


= X 


(18) 


As 


As 


Se 


where  X^^  is  the  mole  fraction  of  As  in  the  liquid  P°g  and  P®^ 


are  the  vapor  pressures  of  As  and  Se  and  P^g  and  Pg^  are  the 
vapor  pressures  in  the  alloy.  Since  the  As  and  Se  vapor  pressures 
are  so  similar  and  the  As  is  only  a small  fraction  of  the  alloy 
composition  this  correction  can  be  neglected. 


C.  Evaporation  Rate 


The  evaporation  rate  in  vacuum  can  be  determined  from  the 
vapor  pressure  using  an  argument  originally  due  to  Langmuir  (see 
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Eq.  32  and  References  29  and  30  for  a critical  discussion) . The 

2 

evaporation  flux  in  g-atoras/cm  sec  is  given  by: 


1 s °v  ^ ^atm 
''  (2nMRT) 


where  M is  the  molecular  weight  in  grams,  MR  = 1.987  cal/g-atom  r 

1/2  1/2  2 

is  the  gas  constant,  f = 1.554  g ' cal  ' /cm  atm  sec  is  a con- 
version factor  and  is  a sticking  or  condensation  parameter. 

This  is  essentially  the  probability  of  evaporation  or  the  proba- 
bility of  condensation.  For  liquid  to  vapor  transitions  a ^ is 
normally  3 1.  The  derivation  of  this  equation  essentially  assumes 

no  interference  from  the  evaporated  material  still  in  the  vicinity 
of  the  evaporating  region.  That  is,  the  derivation  assumes  the 
vapor  pressure  is  much  larger  than  the  background  pressure. 

This  may  not  be  true  at  very  high  evaporation  rates  so  that  Eq.  19 
may  overestimate  the  evaporation  rate  at  high  rates. 


As  an  example,  Eq.  19  gives  for  Se  at  T = 690  °C  where  = 

1 atm,  = 2.64  cm/sec.  Hence  to  machine  0.1  microns  of  material 
at  this  temperature  would  require  3.8  Msec.  To  machine  significant 
amounts  of  material  in-  O.lU  sec  will  require  P~10  atm  or  T-920 
°C. 


. Heat  of  Vaporization 


The  heats  of  transformation  of  Se  (AHj^  for  the  solid  to 
liquid  change  and  Ah^  for  the  liquid  to  vapor  change)  are  tabu- 
lated in  Table  3. 
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TABLE  3 


HEATS  OF  TRANSFORMATION  OF  Se 


Phase  Change 

Reference 

AH^ 

1.25 

kcal/g  atoms 

23 

1.3 

31 

AH 

14.27 

23 

6.3 

31 

It  can  also  be  shown  [ 33  ] that  the  vapor  pressure  is  related 
to  the  heat  of  vaporization  or  the  Gibbs  free  energy  change 
between  the  liquid  and  vapor  state  by: 


log  P 


atm 


4.575T 


+ 


const 


(20) 


Using  the  fit  to  the  vapor  pressure  data  in  Eg.  17  gives  = 

24.6  kcal/g-atoro  which  does  not  agree  well  with  Table  3.  Since 
the  properties  of  high  temperature  Se  have  not  been  well  studied, 
it  is  helpful  to  keep  in  mind  that  the  heats  of  transformation 
in  Table  3 may  be  in  error.  For  the  purposes  of  modeling  we  will 
use  1.25  and  14.3  kcal/g  at  for  the  AH  and  ah^ 


To  estimate  the  effect  of  these  heats  of  transformation  on 
the  writing  process,  consider  the  relative  powers  involved  in 
vaporizing  a (O.lM  volume  of  Se  in  0.1  psec.  Using  the  bulk 

3 

Se  mass  density,  AH^  « 14.3  kcal/g-atom  » 3250  J/cm  . To  vaporize 
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the  O.lp  cube  of  Se  requires  a volume/time  * 3.25  x lO”^  watts 

into  heat  of  vaporization.  As  will  be  seen,  to  maintain  the  tern- 

o ^4 

perature  of  this  volume  of  material  at  900  C requires-  5 x 10 
watts  so  the  effect  of  evaporative  cooling  is  about  a 10%  effect. 
However,  to  machine  the  same  amount  of  material  in  0.01  usee 
would  require  10  times  the  power  into  heat  of  vaporization  which 
is  comparable  to  the  power  required  to  maintain  the  temperature. 
Hence  the  beam  power  would  have  to  be  almost  doubled  to  provide 
the  heat  of  vaporization.  The  heat  of  vaporization  will  become 
a critical  limitation  on  the  writing  rate  above  10  Mbits/sec. 

E.  Backscattered  Electrons 


Se  is  a material  of  relatively  high  atomic  number  so  that 
the  reflected  energy  due  to  backscattered  electrons  is  non- 
negligible.  Figure  35  shows  the  number  backscatter  ration  versus 
atomic  number  Z from  Bishop  [26]  and  Palluel  [27  ] at  5 keV. 

The  energy  dependence  above  5 keV  is  very  weak.  The  energy  back- 
scatter ratio  f is  the  ratio  of  the  total  energy  backscattered 
from  the  sample  material  to  the  incident  electron  energy.  The 
average  energy  of  the  backscattered  electrons  is-  0.5.  A more 
exact  relation  is  given  by  Sternglass  (28]  based  on  integrals 
of  the  backscatter  energy  distributions  as: 

f = (0.45  + 0.002Z)  n (21) 

From  Figure  35  f(Se)  * 0.20  using  Bishop’s  value  for  n • 
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Figure  35.  Number  Backscatter  Ratio  for  5 keV  Electrons 
Versus  Atomic  Number  Based  on  the  Results 
of  Bishop  and  Palluel 


F.  Heating  Models 


As  described  in  a separate  section  of  this  report  a full 
two-dimensional  heating  model  is  being  applied  to  the  calculation 
of  the  temperature  rises  in  multilayer  targets.  In  this  section 
approximate  models  for  the  heating  and  for  the  process  using 
one-  and  two-dimensional  heating  models  will  be  described.  These 
are  combined  with  the  vaporization  properties  in  the  following 
section  to  provide  a complete  milling  model. 
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Figure  36.  One-Dimensional  Heating  Model 
Geometry 


i) One-dimensional  Model 

In  cases  where  the  thickness  of  the  films  is  small  compared 
to  the  electron  range  and/or  electron  spot  size,  one-dimensional 
models  will  provide  very  good  solutions  for  the  temperature  rise. 
Figure  36  shows  the  definition  of  symbols  for  a three  layer  heat- 
ing model.  The  source  of  heat  P(Z)  (electron  beam)  is  uniform 
in  r but  can  vary  arbitrarily  in  depth  z.  This  energy  distribu- 
tion will  be  the  Grun  function  in  the  final  results.  Film  #1 
will  normally  be  Se  and  film  #2  a low  mass  density  thermal  in- 
sulating film  such  as  Si02>  Film  13  will  be  the  silicon  sub- 
strate which  is  taken  as  an  infinite  conductivity  material  for 
this  model  since  its  conductivity  is  about  100  times  Se  or  SiO~. 


- , ....  V 


0 


and  K2  are  the  thermal  conductivities  in  the  two  top  layers. 
The  time  independent  heat  equation  in  1-d  is: 


— - T(Z)  * - P(Z)  - - P,  6(Z) 
K *^1  ^ 


(22) 


where  K is  the  thermal  conductivity  which  is  a function  of  depth 
and  T(Z)  is  the  temperature  distribution.  Pj^  is  the  heat  lost 
from  the  top  surface.  Pj^  will  be  modeled  as  a vapor  cooling 
effect,  although  for  higher  temperatures  a radiation  cooling 
term  should  also  be  included.  This  problem  can  be  solved  by 
Green's  function  techniques,  in  which  the  solutions  in  each  layer 
are  linear.  The  final  solution  for  the  temperature  at  the  top 
surface  is: 

r ^ 


T(o)  = / dZ^  P(Zq) 
o 


^l"*o  . ^2“^1 
*»  • — » 


2 O 


•'■c-e- 


(23) 


ii) Two-dimensional  Solution 


A second  solution  is  for  the  case  shown  in  Figure  37.  This 
geometry  assumes  a semi-infinite  medium  where  the  input  power 
is  uniformly  distributed  over  a hemisphere  of  radius  a.  This 
is  an  appropriate  approximation  when  the  film  thickness  is  much 
larger  or  comparable  to  the  electron  range  and  spot  size.  A 
complete  time  dependent  solution  for  this  problem  was  obtained 
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Figure  37.  Two-Dimensional  Heat- 
ing Model  Geometry 


by  Goldenbecg  [34]  . The  equilibrium  temperature  rise  at  the 
surface  of  the  center  of  the  hemisphere  T(0)  is  given  by: 

T{0)  = 3 P^/4„Ka  (24) 

where  P,j,  is  the  total  power  into  the  hemisphere.  The  thermal 

time  constant  is  given  by; 

I 

, T=  a^k  (25) 

where  k = K/C  is  the  thermal  conductivity.  For  Se  k~4  x 10  ^ 

P 

I cm^/sec  so  that  typical  thermal  time  constants  for  0.1  wm  geometries 

I are  -6  x 10~®  sec.  Hence  we  are  justified  in  treating  all  calcula- 

tions in  the  time  independent  limit  down  to  about  3 thermal  time 

— fl 

constants  of  -2  x 10  sec.  This  was  already  implicitly  assumed 
in  the  one-dimensional  model  discussion. 
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iii)  Effect  of  Electron  Scattering 


The  effect  of  electron  scattering  and  spot  size  will  be 
handled  in  a slightly  different  fashion  for  the  one-  and  two- 
dimensional  models.  For  the  two-dimensional  model,  the  electron 
beam  energy  will  be  distributed  uniformly  in  a disk  of  diameter 
2/3Rg  + b where  b is  the  electron  spot  size  (defined  here  as 
the  1/e  spot  size  i.e.,  the  diameter  where  the  Gaussian  beam 
profile  has  fallen  to  1/e  of  its  peak  value  which  corresponds 
to  the  spot  which  includes  60%  of  the  beam  current) . This  is 
known  to  be  a good  assumption  and  was  discussed  briefly  in  the 
Phase  I First  Quarterly  Report  (April  15  - July  15,  1976)  and 
more  completely  in  Reference  35.  Hence  the  input  beam  power 
distribution  is  given  by; 


(l-P)E„Io 

= -A  R - --  x(Z/Rg) 

Cj 


(26) 


where  Eg  is  the  beam  energy,  Ig  the  total  beam  current,  f the 

energy  backscatter  fraction  discussed  earlier,  x (q)  the  normal- 

2 

ized  Grun  function  Eq.  16  and  A = ti/4  (2/3  Rg+b)  is  the  area 
of  the  film  heated  by  the  beam. 


Equation  26  can  be  generalized  for  the  multilayer  problem  as 
follows; 

r 


P(Z)  = 


X (Z/R^) 

^ ‘(Jr*?) 


V AR2 


Z < Zi 


Z > Zi 


(27) 
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where  and  R2  are  the  Griin  ranges  in  films  1 and  2 respectively 
This  is  a reasonable  assumption  that  was  previously  discussed 
in  connection  with  the  determination  of  the  milled  Se  thickness 
using  electron  beam  diode  signals. 

For  the  two-dimensional  model  the  input  beam  power  is 
given  by: 


= (l-f)Egl3  (28) 

and  the  radius  a of  the  heated  area  by: 

a^  = I (2/3RQ+b)  Rq  (29) 

Equation  29  is  the  geometric  average  of  the  depth  of  beam  penetra 
tion  R_  and  the  lateral  radius  of  the  heated  volume  (2/3R_+b)/2. 
These  two  temperature  models  define  the  extreme  cases  for  the 
heating  model.  Typically  the  one-dimensional  model  conditions 
will  be  well  satisfied  near  the  bottom  of  a milled  hole.  As 
a criteria  for  the  transition  between  the  models  we  require  that 
the  temperature  at  the  surface,  as  a function  of  film  thickness 
Z^,  be  continuous.  This  typically  produces  a transition  between 
the  two  models  for  2/2Ry.  This  is  a reasonable  way  to  join 
the  two  models  that  produces  a slightly  high  estimate  of  the 
surface  temperature  in  the  transition  region.  Hence  the  milling 
rates  will  be  slightly  optimistic. 


iv)  Milling  Model 


In  order  to  treat  the  milling  process  new  depth  variables 
are  defined  as  follows:  Z is  the  depth  of  milled  film  #1  relative 
to  the  initial  thickness.  Let  and  Y2  be  the  thicknesses  of 
the  initial  layers  1 and  2.  Then  = Y^^-Z  and  z^  - Y2+Yj^-Z. 


Using  Eguations  27  and  23  the  one-dimensional  model  for 
the  surface  temperature  can  be  expressed  as 


f-Rj  . z -z  \ 


_p  h + 

*^1  Kj^  Kj 


(30) 


G--  ■''<*5’]  * ^ [H(Zg)  - H(Z5)] 


a Rj>  *1 


G = 0 


aR^i  Zi 


where  H and  F are  particular  integrals  of  A : 


^ 4 q i 

H(q)  = / dyA/y  = L j 

o i = l 


(31) 


" / ydy  A (y)  = z ai_2  ? 


(32) 


1*  1.1  is  the  maximum  extent  of  the  normalized  Grun  function x and 


Zg  » Min  (Zj^,  a Rj^) 
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A 


(33) 


(34) 


r 


Zg  = Min(Z2,  a R2  + Zi  (I-R2/R1)) 

The  term  is  the  power  density  lost  to  heat  of  vaporization. 
This  can  be  expressed  as: 


(35) 


where  is  the  evaporation  flux  in  cm/sec  of  solid  material  and 
A the  heat  of  vaporization  joules/g-atom,  and  and  the 
mass  density  and  g-atomic  weight  of  the  evaporating  layer.  The 
diameter  of  the  evaporating  region  is  effectively  obscured  in 
this  treatment  by  the  way  the  surface  temperature  is  calculated. 
Actually,  only  the  milling  of  the  center  of  the  heated  region 
is  treated  correctly.  Equation  35  assumes  that  the  entire  heated 
area  of  diameter  ~ (2/3Rj^  +b)  is  evaporating  at  the  same  rate. 
Actually  the  edges  will  evaporate  somewhat  slower  so  this  estimate 
of  the  heat  lost  to  vaporization  is  somewhat  too  large. 

Since  Pj^  in  Eg.  35  depends  on  the  surface  temperature  T^ 
through  the  dependence  of  on  T^,  T^  must  be  determined  by 
iteration.  In  the  computer  program  Eqs.  35  and  30  are  solved 
for  T^  by  Gaussian  iteration.  The  temperature  converges  to  with- 
in 0.3%  in  three  iterations. 


The  milling  depth  Z as  a function  of  time  is  determined 
by  integration  of  the  milling  rate  equation  36 


Z(t) 


= ;j^(Vdt 


(36) 
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The  time  step  size  A t for  the  integration  is  adjusted  so  that 
the  incremental  milling  depth  AZ  per  time  increment  is  never 
larger  than  0.003  of  the  first  film  thickness.  This  was  verified 
to  give  good  convergence  of  the  integration. 

RESULTS 


Figure  38  shows  a typical  output  of  the  calculation.  For 
the  moment  ignore  the  experimental  points  and  the  dotted  curves. 

The  vertical  axis  is  the  amount  of  milled  depth  versus  writing 
time  per  spot.  Figure  39  shows  the  surface  temperature  versus 
milled  depth.  This  is  not  the  temperature  profile  in  the  film 
but  a plot  of  the  surface  temperature  as  a function  of  milling 
depth.  The  region  of  constant  temperature  is  the  region  in  which 
the  two-dimensional  model  is  applied.  The  top  curve  is  the  surface 
temperature  without  the  effect  of  heat  of  vaporization  cooling 
and  the  lower  curve  is  the  actual  temperature  during  milling. 

Note  that  the  moderate  effect  of  vapor  cooling  is  ~ 20%  as  pre- 
viously estimated. 

The  solid  curve  in  Figure  38  is  for  100  nA  beam  current 
and  corresponds  approximately  to  the  experiment  described  in 
the  experimental  results  of  Table  1.  The  two  dotted  curves  are 
similar  calculations  for  50  and  200  nA  beam  currents.  This  shows 
the  very  strong  dependence  of  the  writing  rate  on  beam  current 
due  to  the  strong  dependence  of  vapor  pressure  on  temperature. 

Three  sets  of  experimental  points  are  shown  for  140,  70,  and 
50  nA  beam  current.  The  140  and  70  nA  data  comes  from  Table  1 
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while  the  50  nA  data  is  based  on  estimates  of  the  writing  depths 
in  the  secondary  emission  micrograph  in  Figure  40.  In  all  cases 
the  beam  voltage  was  16  keV  and  well  focused  to-  0.1  ^m.  Because 
of  the  large  amount  of  electron  scattering  at  16  keV  the  exact 
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Figure  38.  Milling  Depth  Versus  Time  Calculated  Using 
Milling  Model 

spot  size  is  not  important  since  electron  scattering  determines 
the  heated  volume.  Note  that  for  writing  times  less  than  10~^sec 
there  is  excellent  agreement  between  model  and  experiment.  The 
two  regions  corresponding  to  the  two  temperature  models  can  be 
clearly  distinguished.  The  solid  and  cross  points  for  long  dwell 
times  are  in  the  region  of  the  one-dimensional  model  where  the 
temperature  rise  is  limited  by  heat  flow  to  the  substrate,  while 


Heat  of  UaporizatioA  ■ 0 


Eltctron  bean:  1.ME-M7  Anps  lfi.00  kcU  0.100  nicroits  di«n  l<'c  siMt 
Sclcnian  filn  thicknesa  1.000  nicrons  oxide  thickness  0.100  nicrons 
Scicnittfi  dcnsitji*  4.30s/cnt3  ranse*  i.27nicrons 
th.  cond.-O.OOSOwatts/’cntZ  de*  heot  vap.«  5.97E+004  j/f-ot 
Initial  Surface  Tenp  rel  to  anbient*  766  desrecs 

Figure  39.  Surface  Temperature  Versus  Residual  Se  Depth  for 
the  Case  in  Figure  38.  Top  curve  is  without 
the  heat  of  vaporization  correction. 


at  the  same  time  the  open  circle  points  are  in  the  infinite  media 
approximation  region  where  the  temperature  surface  temperature 
is  independent  of  film  thickness  and  is  determined  by  the  thermal 
conductivity  of  the  Se  film. 


As  can  be  seen,  the  agreement  is  very  good.  A further  improve- 
ment could  be  obtained  by  multiplying  all  of  the  calculated  surface 
temperatures  by  about  1.3  since  in  all  cases  the  data  is  better  fit 
by  a beam  current  about  1.3  times  larger  than  measured.  The  small 
difference  between  theory  and  experiment  is  either  due  to  the  ap- 
proximate thermal  models  used  or  to  a lower  thermal  conductivity 
within  the  range  of  published  values.  Since  the  error  is  such  as 
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SEM  Secondary  Image  of 
Writing  on  1.8y  Se  at 
50  nA  and  16  keV  for  In 
dicated  Dwell  Times 


to  overestimate  the  required  beam  current  no 
or  refinements  will  be  attempted  at  this  time 


Note  that  the  model  deviates  from  experiment  at  long  dwell 
times.  This  is  due  to  melting  and  flow  of  the  film.  The  melting 
temperature  is  200  °C  which  is  a factor  of  3 below  the  temperature 
required  for  rapid  vaporization.  Observations  of  the  Se  film 
in  the  SEM  under  raster  scan  conditions  show  that  as  the  Se  film 
is  raised  above  its  melting  temperature  by  the  beam  it  tends 
tc  wet  and  flow  towards  the  region  last  hit  by  the  electron  beam. 
Hence,  if  the  beam  dwells  for  a long  time  on  one  spot  molten 
material  at  the  bottom  of  the  milled  hole  will  flow  up  the  side 


of  the  walls  increasing  the  milled  depth.  The  observed  flow 
rates  are  very  slow.  This  is  consistent  with  the  very  high  vis- 
cosity of  molten  Se.  This  explanation  has  been  confirmed  by 
writing  with  much  lower  beam  currents  with  power  sufficient  only 
to  melt  the  top  surface.  Melting  and  flow  of  the  film  to  the 
sides  of  the  heated  region  in  -1  msec  were  observed. 

With  this  justification  it  is  now  appropriate  to  apply  the 
model  to  predictions  of  the  required  beam  current,  spot  size, 
and  voltage  for  high  density  writing.  To  achieve  0.1  ym-bit 
sizes,  it  will  be  necessary  to  produce  heated  volumes  of  dimen- 
sion 0.1  ym.  Because  of  electron  beam  lateral  scattering,  this 
will  require  beam  voltages  less  than  5 keV  and  spot  sizes  slightly 
less  than  0.1  ym.  For  this  reason  the  model  calculations  described 
in  the  remainder  of  this  section  concentrate  on  parameters  in 
this  region.  Figures  41  through  45  show  some  applications  of 
the  model  to  such  high  density  writing  conditions.  All  of  the 
relevant  parameters  can  be  read  in  the  headings.  The  total  diam- 
eter of  the  heated  region  will  be  approximately  the  spot  diameter 
+ 0.66  times  the  electron  range  in  Se.  The  election  range  in 
Se  is  indicated  in  the  heading. 

Figures  41  and  42  show  the  milling  behavior  for  a 3 keV 
electron  beam  on  a 0.03  ym  thick  Se  layer  with  0.05  ym  of  oxide 
for  thermal  isolation.  A beam  current  of  100  nA  removes  0.024  ym 
of  Se.  The  residual  amount  of  Se  is  0.06  ym.  If  a 3 keV  electron 
beam  is  used  for  readout  of  this  milled  region  the  diode  signal 
gain  will  be  approximately  250.  This  is  0.3  of  the  maximum  possible 
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Electro*  been:  5.0OE-088  Anps  3.96  kcU  0.039  nicrons  dien  l<'e  siMit 
Scicoittii  filn  thickness  9.039  nicrons  oxide  thickness  9.999  Microns 
SeleniuM  densitu^  4.399''CMt3  ransc*  9.97fiicrons 
th.  cond.>9.99S9watts/cnt2  des  heat  vap.«  S.97E^994  j/f'et 
Initial  Surface  Tenp  rel  to  unbient*  1996  desrees 

Figure  41.  Milling  Calculation  for  3 keV,  50  nA  Q.03U  Beam 
on  0.03u  Se  on  0.05u  Oxide 


gain  that  would  occur  if  there  were  no  residual  Se  and  no  oxide 
isolation  layer.  Since  the  diode  signal  in  the  unwritten  area 


would  be  zero  the  gain  of  250  should  provide  sufficient  signal 
for  10  Mhz  readout  rates  with  beam  currents-  10  nA. 


Figure  43  shows  a similar  calculation  without  the  oxide 
isolation  layer  but  with  the  Se  film  thickness  increased  from 
0.03  to  0.055  uin  to  provide  an  equivalent  mass  thickness  of  ma- 
terial above  the  diode.  At  10~^  sec, 0.028  y m of  Se  is  milled 
as  compared  to  0.024  uni  for  the  above  case  with  the  oxide  isola- 
tion layer  with  the  result  that  the  diode  signal  is  10%  higher. 
Hence  the  oxide  thermal  isolation  layer  has  no  clear  advantage. 


77 


fl.' 


Electro*  bom:  1.ME-M7  Anps  3.M  keV  t.U*  nicrow  4i«n  l^o  «Mt 
So1««i«N  filn  thickness  •.•9t  nicrons  oxitic  thickness  •.•St  elcrene 
Scicninfi  density*  4.399/'ci«t3  mnee*  •.•Teicrons 
th.  cond.*t.ttSt»Mtts/cnt2  dcs  best  vnp.*  S.97Eittt  J/e-nt 
Initial  Sarface  Tenp  ral  to  anbient*  1431  daeraes 

Figure  42.  Milling  Calculation  for  3 keV,  100  nA 
0.03m  Beam  on  0.03ij  Se  on  O.OSy  Oxide 

Of  course  significant  changes  in  the  thermal  conductivities  of 
either  the  milled  or  isolation  layer  would  change  this  conclusion. 


Figure  44  shows  a similar  calculation  at  4 keV  for  the  case 
of  no  thermal  isolation.  The  Se  thickness  has  been  increased 
to  0.11  urn  which  is  the  range  of  the  electron  beam  so  that  the 
diode  readout  signal  at  4 keV  will  still  be  approximately  zero. 


Again  100  nA  produces  significant  milling.  Figure  45  shows 
the  surface  temperature  profile  for  the  case  of  Figure  44.  Note 
the  significant  effect  of  the  heat  of  vaporization  that  reduces 
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Electro*  taoN:  1.ME-M7  Anps  3.M  kcM  •.•3t  nicrons  4i«ii  1/c  s»ot 
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Figure  43.  Milling  Calculation  for  3 keV,  100  nA 
on  O.OSSp  Se 


the  surface  temperature  by  almost  a factor  of  two.  This  is  typ- 
ical of  writing  cases  where  the  surface  temperature  is  such  as 
to  produce  significant  milling  in  less  than  0.1  ysec.  Figure  46 


illustrates  the  effect  of  heat  of  vaporization  on  writing  rate. 

In  the  absence  of  heat  of  vaporization,  doubling  the  beam  current 
should  increase  the  writing  rate  by  several  orders  of  magnitude 
because  of  the  strong  dependence  of  evaporation  rate  on  tempera- 
ture. This  phenomenon  is  seen  to  occur  for  beam  currents  between 
50  and  100  nA.  However,  at  the  higher  currents  almost  no  increase 
in  writing  rate  is  obtained  because  most  of  the  added  power  is 
going  into  heat  of  vaporization.  Hence,  the  large  heat  of  vapor- 
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Figure  44.  Milling  Calculation  for  4 keV,  100  nA 
0.05 u Beam  on  O.lly  Se 


ization  of  Se  will  effectively  limit  the  writing  rate  to  about 
50  MHz  for  beam  currents  300  nA. 


Figure  47  shows  the  milling  behavior  predicted  for  a 25  nA, 
0.1  ijm  diameter  beam  at  5 keV.  This  is  comparable  to  the  best 
performance  currently  available  at  low  voltage  in  the  C&W  writing 
station.  Occasionally  higher  currents  are  achieved  but  at  in- 
creased spot  sizes  so  that  the  beam  power  density  is  actually 
reduced.  No  milling  is  predicted  for  dwell  times  below  10~^  sec. 
However,  as  previously  discussed,  viscous  flow  of  the  molten 
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Eltctron  bean;  1.0DE-O87  Anps  4.M  kcU  0.996  nicroits  dian  spot 
ScIcsisR  film  thickness  0.110  nicrons  oxids  thickness  0.000  nicrons 
Seleniiifi  density^  4.30p/cnt3  ranpe-  D.llnicrons 
th.  cond.B9.9959watts/'ciit2  dep  heat  vap."  9.97E*004  J<'f-at 
Initial  Surface  Tenp  rcl  to  anbient*  1267  deprees 

Figure  45.  Surface  Temperature  During  Milling 
for  Case  in  Figure  44. 


Se  occurs  for  dwell  times  above  10  sec  making  interpretation 
of  long  dwell  time  experiments  difficult. 

Note  that  the  predicted  surface  temperature  is  just  above 
the  melting  temperature  of  200  °C  for  Se.  Viscous  flow  writing 
is  undoubtedly  the  origin  of  the  writing  shown  in  Figure  31. 

RELATION  BETWEEN  MASS  DENSITY  AND  THERMAL  CONDUCTIVITY 


FOR  THE  CHOICE  OF  AN  OPTIMUM  MILLING  MEDIA 


Low  thermal  conductivity  is  obviously  desired  for  thermal 
milling  because  the  maximum  temperature  rise  is  inversely  proper- 


Electro*  been:  1.ME-M7  Anps  4.M  kcU  9.1M  nicroM  dien  1/c  S90t 
SclcuiiM  filn  thicimess  9.  lit  nicrons  oxide  thickoes*  9*999  nicroM 
Sclcnien  density*  4.399<'cnt3  ranyc*  9. llnicrons 
th.  cond.^t.ttStwatts/cntZ  dey  hc«t  vap.*  5.97E‘»994  j<«'9-ot 
Initial  Surface  Tenp  rcl  to  anbient*  1996  deyrees 

Figure  46.  Milling  Calculation  for  4 keV,  100  nA 
O.Olu  Beam  on  O.lly  Se 


tional  to  the  thermal  conductivity,  ignoring  substrate  effects. 


However,  to  some  extent,  a higher  mass  density  can  compensate 


for  a higher  thermal  conductivity  because  the  power  absorption 
per  unit  depth  from  the  electron  beam  is  greater.  Higher  electron 
beam  energies  can  then  be  utilized  because  of  the  reduced  elec- 
tron range  and  lateral  scattering  in  the  higher  mass  density 
materials.  Based  on  the  two-dimensional  heating  model  discussed 


previously,  (see  Equation  24)  the  surface  temperature  rise  in 
the  infinite  media  approximation  is  given  by: 


i 


Electron  kc«fi:  SeSSE-MS  Anps  5.M  kcU  •. IM  nicrons  dicn  1/e  sect 
Seleeien  file  ttiickness  0.8M  Microns  oxide  thickness  t. IM  Microns 
Seleninn  dcnsiti^  4.309/CMt3  ransc*  B. 17nicrons 
th.  cond.«8.005Bw«tts/cnt2  des  hent  v«p.*  S.BTE^BBd  J/9~et 
Initial  Surface  Tcmp  rcl  to  aMbicnt*  209  desrccs 

Figure  47.  Milling  Calculation  for  5 keV,  25  nA, 

O.lM  Beam 


where  = Ig  x Eg  is  the  beam  power.  This  ignores  spot  size 
effects  so  that  the  radius  of  the  heated  volume  is  taken  as  the 
Grlin  range  Rq.  Using  Eg.  15  for  results  in  Eg.  38  for  the 
relation  between  surface  temperature  at  the  start  of  milling, 
mass  density,  and  thermal  conductivity.  Requiring  a surface 
temperature  of  2000  °C  for  rapid  thermal  milling  and  R^,  = O.lu  to 
limit  the  size  of  the  heated  and  milled  region  to  O.ly  yields 
the  curves  in  Figure  48.  Any  material  has  a unique  position 
on  this  plot.  The  solid  lines  labeled  with  various  beam  currents 
are  the  locus  of  points  for  which  the  surface  temperature  rise 
will  be  exactly  2000  °C  based  on  Eq.  38  for  R^  = O.lp. 
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FILM  MASS  DENSITY  g/cm» 

Figure  48.  Thermal  Conductivity  Versus  Mass  Density  Plot 
for  Various  Materials  Compared  to  a 
Heating  Model  for  the  Maximum  Tem- 
perature Rise  Achievable  in  Various 
Films  at  0.1  Micron  Written  Spot  Size 
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Any  material  below  this  line  will  reach  a temperature  higher 
than  2000  for  the  indicated  beam  current.  This  graph  illus- 
trates that  the  trade-off  between  mass  density  and  thermal  con- 
ductivity goes,  approximately  as  the  square  root  of  the  mass 
density  and  provides  a convenient  method  for  comparing  the  rela- 
tive writing  efficiencies  of  materials  of  different  mass  densities. 
Plotted  on  the  graph  are  data  points  for  some  common  materials. 

The  polyimide  point  is  typical  of  many  polymers.  This  shows 
very  clearly  that  even  though  many  polymer  films  might  have 
slightly  lower  thermal  conductivity  than  Se-As,  their  reduced 
mass  density  does  not  make  them  any  more  attractive.  Of  course, 
if  the  mass  density  could  be  increased  by  loading  with  heavy 
metals  without  seriously  affecting  the  thermal  conductivity  this 
process  would  be  preferable.  Metal  loading  might  also  reduce 
the  other  problem  with  polymers  of  electrostatic  charging  during 
electron  beam  readout  or  writing. 


SUMMARY 


Electron  milling  by  evaporation  has  been  demonstrated  on 
amorphous  Se-As  alloy  evaporated  films.  A model  for  the  milling 
process  has  been  developed  which  includes  the  effects  of  thermal 
evaporation,  solutions  to  the  heat  flow  equation,  and  evaporative 
cooling.  This  has  been  compared  to  writing  experiments  at  16 
and  5 keV  on  0.05  and  1.8  ym  thick  evaporated  Se-As  films.  Good 
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agreement  was  obtained  between  these  experiments  and  the  model. 
Based  on  the  model  it  is  predicted  that  a 0.05  ym  diameter  beam 
at  4 keV  and  100  nA  or,  0.1  ym  diameter  beam  at  200  nA  will  be 
sufficient  for  10  MHz  milling  at  0.1  ym-bit  spacing.  The  above 
beam  specifications  assume  that  all  but  1/2  of  the  beam  current 
is  contained  in  the  stated  beam  diameter.  The  effect  of  evapo- 
rative cooling  is  small  (20%)  but  does  have  a noticeable  effect 
at  10  MHz  writing  rates.  The  writing  rate  will  be  limited  by 
evaporative  cooling  and  the  thermal  time  constant  to  about  50  MHz 
for  beam  currents  less  than  400  nA.  The  Se-As  films  prepared  for 
these  studies  are  amorphous,  uniform  and  stable,  and  have  the 
thermal  properties  expected  based  on  reported  measurements  on  simi- 
lar bulk  materials. 
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SECTION  V 


TWO-DIMENSIONAL  ELECTRON  BEAM  HEATING  MODEL 

INTRODUCTION 

The  increasing  complexity  and  diversity  of  the  multilayered 
electron  beam  targets  under  examination  have  resulted  in  a greater 
need  for  a more  detailed  and  precise  method  of  modelling  the  elec- 
tron beam  heating.  The  simple  analytical  approaches  can  only  be 
used  with  confidence  to  model  single  thick  layers  or  to  gain  a 
basic  understanding  of  the  temperature  increase  on  beeun  diameter 
and  thermal  conductivity.  Our  interest  in  material  systems  using 
two  and  three  materials  with  widely  differing  densities  and  ther- 
mal properties  has  made  it  impossible  to  predict  with  any  confi- 
dence either  the  distribution  of  electron  beam  energy  loss  or  the 
resulting  temperature  patterns.  The  simple  computer  models  de- 
scribed in  the  first  quarterly  report  have  been  expanded  to  a 
family  of  programs  designed  to  deal  with  the  variety  of  needs* 
each  in  the  most  efficient  possible  way.  In  addition*  a model 
for  electron  beam  heat  loss  has  been  formulated.  We  believe  that 
this  model  closely  approximates  an  actual  beam  even  in  a target 
of  multilayered  materials.  The  model  is  based  on  actual  experi- 
mental studies  and  will  be  described  in  detail  later. 

The  very  small  physical  dimensions  of  the  heated  area  and 
the  very  short  times  required  for  the  completion  of  these  processes 
conspire  to  make  it  virtually  impossible  to  examine  any  sort  of 
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details  of  the  process  on  an  exper iemental  level.  Thus  experi- 
mentally^  the  process  either  works  or  it  doesn't.  The  existence 
of  a detailed,  time  dependent  model  allows  us  to  understand  in 
detail  why  a given  process  either  succeeds  or  fails  and  what  changes 
are  likely  to  improve  it.  The  model  also  permits  us  to  try  experi- 
ments that  are  presently  beyond  our  capability  but  that  may  become 
practical  in  the  future.  Finally,  it  permits  us  to  vary  the  large 
number  of  parameters  involved  in  the  system  and  to  determine  the 
optimum  values.  In  much  of  the  discussion  which  follows,  a par- 
ticular case  will  be  used  as  an  example.  The  case  chosen  is  a 

e o 

600  A film  of  selenium  over  a 400  A film  of  silicon  dioxide  on  ’ 

! i 

a silicon  substrate.  The  electron  beam  is  4 kV  at  1 yA  and  has 
a 500  A radius. 

MODEL  CONCERNS  ^ 


No  model,  no  matter  how  carefully  formulated,  can  exactly 
reproduce  nature.  It  is  important  to  recognize  in  advance  what 
approximations  are  being  made  and  attempt  to  assess  their  impact 
on  the  validity  of  the  results. 

1.  Bulk  values  are  assumed  for  all  of  the  material  parameters 
in  spite  of  the  very  small  dimensions  and  the  possibility 

of  size  effects. 

2.  All  interfaces,  whether  between  differing  materials  or 
between  the  material  and  its  outer  boundaries,  are  assumed 
to  be  without  substance  even  though  we  know  that  there  is 
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a different  physical  and  electronic  structure  near  these  in- 
terfaces. This  is  done  simply  because  we  do  not  know  how 
to  treat  them  properly. 

3.  The  model  is  a finite  difference  formulation  of  the  heat 
flow  equation  with  the  electron  beam  energy  loss  function 
providing  the  heat  input.  As  a part  of  the  solution  the  de- 
vice is  broken  up  into  finite  cell  volumes  and  the  heat  flow 
equation  is  approximated  at  these  points  by  finite  differences 
replacing  the  spatial  derivatives. 

The  effect  of  the  finiteness  of  cell  size  is  measured 
by  reducing  cell  size  and  observing  the  difference.  Figure 
49  shows  the  center  temperature  as  a function  of  depth  for 
our  example  after  0.2  ns.  The  open  circles  of  the  coarser 
mesh  show  that  this  faster  and  more  economical  size  is  com- 
pletely adequate  for  most  purposes. 

4.  The  peripheral  cells  are  held  at  ambient  (starting)  tem- 
perature. This  brings  a finiteness  to  the  problem  and  repre- 
sents a "worst  case"  temperature  rise. 

5.  The  program  has  been  tested  by  comparison  with  the  "uniform 
hemisphere"  problem  which  can  be  done  analytically.  The  re- 
sults are  also  constant  with  varying  time  interval  and  cell 
size . 

6.  The  distribution  of  energy  loss  from  the  electron  beam 
has  been  studied  Cor  a number  of  materials  and  we  have  a high 
degree  of  confidence  with  the  present  algorithm.  One  lingering 
concern  is  that  energy  is  absorbed  by  a semiconductor  pri- 
marily by  electron-hole  excitation  and  that  the  actual  phonon 
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(heat)  production  may  occur  at  a different  location  depending 
on  drift  fields,  lifetime,  etc. 


MODEL  DESCRIPTION 


f I 
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The  most  general  version  of  the  computer  program  has  [1] 
variable  cell  size  [2]  checks  every  cell  boundary  for  continuity 
of  heat  flow  in  the  event  of  differing  adjacent  materials  i3J, 
and  takes  account  of  the  temperature  dependence  of  all  thermal 
properties.  The  most  useful  version  has  constant  cell  size.  That 
is,  the  cross-sectional  area  of  the  ring-shaped  cells  is  fixed. 

The  cell  volume  obviously  depends  on  the  "ring"  diameter,  and  the 
radial  (R)  and  axial  (Z)  dimensions  need  not  be  the  same.  This 
version  assumes,  at  most,  two  uniform  layers  of  different  material 
on  a still  different  substrate.  It  also  suppresses  the  tempera- 
ture dependence  of  the  thermal  characteristics  (including  phase 
changes)  in  order  to  both  speed  the  calculation  and  permit  direct 
— temperature  scaling  with  beam  current.  The  following,"  then,  is 
a description  of  this  most  useful  program,  "H0TSP0T3."  A short 
user's  guide  is  presented  in  Appendix  A followed  by  a listing  of 
the  program. 

The  development  of  a single  general  purpose  algorithm  to  de- 
scribe the  electron  beam  heating  has  required  considerable  effort 
and  trial.  The  first  step  was  to  solve  the  problem  in  a single 
medium.  Fortunately,  an  empirical  formula  for  total  heat  loss 
as  a function  of  depth  has  already  been  developed  and  depends 
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Figure  49.  Axial  Temperature  as  a Function  of  Depth 
Comparing  Two  Different  Mesh  Sizes 


only  on  density  of  the  target  material  and  beam  energy.  A uni- 
versal energy  loss  per  unit  volume  curve  can  be  drawn  if  the 
depth  argument  is  given  in  distance  divided  by  Griih  range,  Rg. 
"In  order  to  extend' this  "universal  curve'tr(y/Rg)  “from  one  material 
to  another,  it  was  assumed  that  the  normalized  depth  argument, 
Y/Rg,  was  continuous  at  boundaries  between  materials  of  differing 
range.  The  radial  shape  of  the  beam  is  arbitrarily  taken  to 
be  Gaussian  at  all  depths.  The  only  remaining  problem,  then, 
is  to  decide  in  what  way  the  width  of  the  Gaussian  varies  with 
depth.  After  many  trials,  it  was  found  that  best  agreement  with 
the  experimental  data  available  was  obtained  if  the  width  in- 
creased linearly  at  a 30  degree  angle  from  the  beam  axis. 
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Figure  50  shows  a plot  of  L(y/Rg)  as  a function  of  depth  in 


our  multilayered  example.  Also  shown  is  Y/Rg  (dashed  line).  Fig- 
ures 51  and  52  show  contours  of  equal  energy  loss  per  volume  for 


beams  differing  only  in  energy.  These  diagrams  illustrate  one 


' of  the  subtleties  of  this  problem  and  further  demonstrate  the 
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need  for  such  a model.  One  might  expect  that  simply  increasing 


the  beam  energy  would  increase  all  of  the  temperatures  in  the  spot 


the  increase  in  volume  of  the  beam  pattern  with 


area.  However 


increasing  energy  can  more  than  offset  the  greater  available  power 
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BEAM  SHAPE  Te 
7.0  KeV 


Contours  of  Equal  Energy  Loss  for  a 1000  A Diameter 
7 kV  Electron  Beam  in  Tellurium 


Figure  52 


In  fact,  there  is  an  optimum  beam  energy  for  maximizing  the  tem- 
perature at  any  given  depth.  Such  matters  will  become  clearer 
during  the  discussion  of  results. 


SELECTED  RESULTS 


Some  examples  will  now  be  given  of  uses  to  which  the  program 
has  been  put.  For  purposes  of  scaling  convenience,  a beam  cur- 
rent of  one  microampere  has  been  used  throughout  even  though 
this  leads  to  unrealistic  temperatures  in  some  cases. 

1.  It  was  desired  to  learn  the  effect  of  a thin  layer  of 
silicon  dioxide  between  the  silicon  substrate  and  a 600  A 
selenium  film.  This  could  be  useful  for  electron  beam  milling 

O 

The  program  was  run  with  and  without  a 400  A silicon  dioxide 
middle  layer.  Figure  53  shows  the  axial  temperature  as  a 
function  of  depth  after  1.5  ns.  The  relatively  low  thermal 
conductivity  of  the  oxide  layer  has  resulted  in  a somewhat 
higher  temperature  in  the  selenium.  If  the  program  is  ex- 
tended to  4 ns  it  can  be  seen  from  Figure  54  that  the  effect 
of  the  oxide  is  relatively  much  greater  at  longer  heating 
times. 

Another  concern  is  the  size  of  the  mark  made  on  the  sur- 
face. To  examine  this,  a family  of  curves  showing  the  tempera 
ture  at  various  depths  as  a function  of  distance  from  the 
beam  axis  is  shown  in  Figure  55.  Note  that  the  diameter 
affected  may  be  either  larger  or  smaller  than  the  effective 


beam  diameter  of  500  A depending  on  where  the  vaporization 
temperature  occurs  on  each  curve.  In  all  cases  the  beam 
diameter  is  defined  as  the  e^  diameter  of  a Gaussian  spot 
profile. 

2.  In  some  larger  scale  milling  experiments,  it  was  decided 
to  calculate  the  surface  heating  rate  of  a 12000  X selenium 
film  on  silicon  dioxide  using  a 16  kV  beam.  Figure  56  shows 

O 

the  surface  (500  A deep)  temperature  as  a function  of  time 
plotted  for  two  different  beam  radii. 

3.  In  a study  of  the  effect  of  beam  energy  and  film  thickness 
on  the  temperature  of  tellurium  films  on  silicon,  the  follow- 
ing results  were  obtained  (see  Table  4) . Over  the  range 

O 

considered,  the  interface  temperature  for  a 600  A tellurium 
film  on  silicon  goes  through  a maximum  at  a beam  energy 
between  5 and  6 kV.  Yet  a still  higher  interface  temperature 

O 

can  be  reached  with  a thinner  (400  A)  film  at  4 kV. 

4.  Calculations  done  for  gold  and  aluminum  on  silicon  have 

o e 

also  been  performed.  A 4 kV  — A diameter  beam  on  a 200  A 
gold  film  will  raise  the  interface  temperature  by  72  °C 

O O 

in  0.2  ns.  A 2.5  kV  -A  diameter  beam  on  a 465  A aluminum 
film  will  raise  its  interface  temperature  by  60  °C  in  0.6  ns. 
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DEPTH  FROM  SURFACE  (ANGSTROMS) 

Axial  Temperature  as  a Function  of  Depth 
After  1.5  ns  Showing  Influence 
of  the  Silicon  Dioxide  Layer 


Figure  53 


TIME  (nano  sk) 


Figure  54.  Surface  Temperature  (100  A deep) 
as  a Function  of  Time  Showing 
Influence  of  the  Silicon  Dioxide 
Layer 
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TEMPERATURES  AT  VARIOUS  DEPTHS 


DISTANCE  FROM  BEAM  AXIS  -ANGSTROMS 


Figure  55.  Temperatures  Away  from  the  Beam  Axis 
at  Various  Depths  After  5 ns 
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Figure  56.  Surface  Temperature  (500  A deep) 
as  a Function  of  Time  Showing 
the  Influence  of  Electron  Beam 
Radius 


TABLE  4 

CALCULATED  TEMPERATURES  AFTER  1.5  ns 

O ® . • 

For  a — A diameter  beam  on  a 600A  Te  film  on  Si 
Beam  Energy  Te  Film  Temp,  deg  K Temp,  deg  K Temp,  deg  K 


kilovolts  angstroms 


100  A 


1169 

1039 


300  A 


interface 


SECTION  VI 


I: 


ELECTRON  MICROPROBE  ANALYSIS 

INTRODUCTION 

During  the  current  phase  of  the  Archival  Memory  Contract, 
considerable  interest  in  electron  beam  thermal  machining  as  a re- 
cording technique  has  developed.  With  this  method  of  recording, 
a high  current  finely-focused  electron  beam  is  used  to  machine 
submicron  data  sites  on  the  surface  of  a suitable  recording  medium. 
Electron  beam  thermal  machining  as  a method  of  recording  informa- 
tion is  not  a new  concept  in  itself  and  has  been  previously  worked 
on  at  GE  as  well  as  in  other  companies.  [ 37/381  Prior  efforts 
in  this  area  tended  to  be  limited  in  target  material  and  operation 
due  to  the  low  current  capabilities  of  submicron  electron  probes. 

A further  limitation  was  the  requirement  for  secondary  emission 
readout.  During  recent  years  the  thermionic-aided  field  emission 
(TFE)  electron  source  has  seen  significant  advances  in  its  develop- 
ment to  the  point  where  now  it  can  be  considered  a reliable  high 
current  density  source  for  submicron  operation.  This  source  coupled 
with  the  development  of  efficient  diode  detectors  in  both  the  BEAMOS 
and  Archival  Memory  Programs  makes  the  thermal  machining  memory 
concept  a much  more  attractive  technology. 

A design  for  a TFE  source  electron  probe  capable  of  one  micro- 
amp in  a 1000  % spot  (90%  current  diameter)  was  developed  during 


Phase  I of  this  contract  and  reported  in  the  Phase  I Fourth  Interim  |* 

Report  (January  16  - April  15,  1977.  Subsequently,  it  was  learned  “ 

that  Lee  Veneklassen  of  ETEC,  one  of  the  leaders  in  TFE  development, 
was  developing  a similar  electron  probe.  Negotiations  began  earlier  | 

in  this  reporting  period  with  ETEC  for  the  purchase  of  a column  j 

of  this  type  for  a research  tool.  In  order  for  us  to  better  under- 
stand the  column  specifications  and  its  projected  capabilities, 
an  analysis  of  the  microprobe  optics  was  performed  under  this  con- 
tract. The  analysis  is  described  in  this  section  of  the  report. 

MICROPROBE  ANALYSIS 

By  way  of  background  the  column  schematic  for  the  electron 
microprobe  design  developed  here  during  Phase  I is  shown  in  Fig-  j 

ure  57.  As  can  be  seen,  the  microprobe  consists  of  a thermionic- 
aided  field  emission  source  with  a dual  magnetic  lens  focusing 
system.  The  design  relied  upon  the  first  lens  collimating  the 
beam  and  the  second  lens  producing  the  final  image  on  the  target. 

In  order  to  obtain  low  aberrations  and  therfore  much  higher  current 
I density,  the  target  plane  had  to  be  close  to  the  final  lens.  There- 

fore, prelens  double  deflection,  as  indicated  in  this  schematic, 

, was  required.  Based  on  the  analysis  of  this  design,  the  resulting 

current  versus  spot  size  characteristics  as  shown  in  Figure  58 
: were  developed.  As  seen  from  this  figure,  at  4000  V,  better  than 

I O 

one  microamp  of  current  can  be  delivered  into  a 1000  A spot.  The 
analytical  techniques  combine  aberrations  and  geometric  spot  sizes 
I i<^  quadrature.  Previous  experience  indicates  that  this  represents 
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a 90%  current  diameter  spot.  The  proposed  ETEC  design  as  published 
earlier  this  year  39  is  shown  in  Figure  59,  along  with  experi- 
mentally obtained  current  versus  spot  size  data.  In  this  column 
the  electron  gun  is  an  oxygen  processed  100  thermionic-aided 
field  emission.  A specially  designed  variable  focus  electrostatic 
lens  serves  as  the  first  lens.  The  final  lens  is  a standard  mag- 
netic ETEC  autoscan  objective  lens  with  prelens  deflection.  For 
the  high  current  density  probe  it  is  seen  that  a slightly  converg- 
ing probe  between  lenses  is  used  rather  than  a collimated  system 
as  shown  in  our  design.  The  current  versus  spot  size  characteris- 
tic, however,  indicates  that  for  the  high  current  mode  about  one 
microamp  of  current  is  achievable  in  the  tenth  micron  spot  diameter. 
In  fact,  the  whole  probe  current  versus  spot  size  characteristic 
appears  very  similar  to  the  Phase  I design  shown  previously  in 
Figure  58.  With  this  data  negotiations  began  with  ETEC  for  a simi- 
lar electron  probe  system. 

In  order  that  we  might  better  understand  the  ETEC  column  de- 
sign and  its  specifications,  an  analysis  of  the  focusing  optics 
was  conducted  in  this  Laboratory.  Optical  properties  for  the  aber- 
ration coefficients,  focal  lengths,  etc.  were  requested  from  ETEC 
for  their  column  design.  These  properties  were  supplied.  Using 
this  data  in  our  analytical  programs,  the  beam  current  versus  spot 
size  characteristic  shown  in  Figure  60  was  calculated.  The  fact 
that  these  computations  resulted  in  a probe  current  of  only  a few 
tens  of  nanoamperes  in  a 1000  A spot  size  was  totally  unexpected 
and  quite  perplexing.  On  the  one  hand,  it  was  recognized  that 
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Figure  57.  Electron  Microprobe  Design  of  Phase  I 

Veneklassen ' s data  was  experimental  and  thus  hard  to  refute.  This 
made  the  calculation  procedures  used  here  of  some  question.  On 
the  other  hand,  similar  calculation  procedures  have  previously 
been  used  in  BEAMOS  and  other  programs  at  General  Electric  labora- 
tories with  excellent  correlation  between  experiment  and  theory. 

The  major  responsibility  of  this  analytical  study  was  then  to  re- 
solve the  discrepancy  between  these  initial  calculations  and  the 
experimental  results  published  for  the  column  similar  to  the  column 
that  was  in  negotiations. 

One  of  the  first  possible  areas  for  disagreement  was  noted 
to  be  in  the  definition  of  spot  size.  In  the  experiments  of 
Veneklassen,  spot  size  was  determined  by  scanning  the  beam  across 
a sharp  knife  edge.  The  resulting  rise  time  and  current  achieved 
by  scanning  across  the  sharp  edge  were  then  correlated  with  the 
scan  velocity  to  result  in  a given  distance  measurement  for  the 
15-85%  rise  time.  As  previously  noted,  the  spot  size  for  the 
analytical  work  was  determined  by  combining  the  various  spot  sizes 
due  to  the  aberrations  and  the  geometric  imaging  of  the  source 
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Figure  59. 


ETEC  - Column  and  Experimental  Beam  Current 
vs.  Spot  Size  Data 


in  quadrature.  The  spot  size  determined  in  this  manner  corresponds 
to  the  diameter  which  contains  90%  of  the  total  probe  current. 
Depending  on  the  distribution  function  for  the  current  density 
of  the  electron  beam  in  the  image  spot,  the  difference  between 
the  90%  containing  spot  diameter  and  the  15-85%  rise  time  measure- 
ment can  be  considerably  different.  To  show  this,  consider  the 
generally  assumed  Gaussian  current  density  distribution  as  given 
by  Eg.  39: 


J(r)  = 


° - (-^)2 
-o  e 'r_  ' 


(39) 


Tf  r. 


where  1^  is  the  total  current  in  the  spot 
r is  the  spot  radius 


r is  the  radius  with  the  beam  at 


The  current  contained  with  a radius  of  spot  r^ 
this  equation  by  integrating  with  the  result: 


can  be  found  from 


(40) 


The  current  scanned  over  a knife  edge  is  the  convolution  of  the 
current  distribution  given  by  Eg.  39  and  the  unit  step  represent- 
ing the  knife  edge.  The  results  of  such  a convolution  are  given 
in  Eqs.  41  and  42: 


I 


(41) 


(42) 


The  knife  edge  and  the  resulting  rise  time  scan  is  shown  graphi- 
cally in  Figure  61.  Calculation  of  the  spot  size  for  various  rise 
times  based  on  the  convolution  Eqs.  41  and  42  and  correlation  with 
the  current  contained  with  a spot  of  radius  80,  90  and  99%  as  de- 
termined from  Eq.  40  are  given  in  Table  5. 


! 


As  can  be  seen  from  the  data  in  Table  5,  the  resulting  90%  cur- 
rent profile  for  Gaussian  beam  is  2 times  the  measured  spot  size 
as  determined  for  15-85%  rise  time.  Therefore,  a factor  of  2 can 
be  accounted  for  due  to  the  difference  in  the  way  the  spot's  size 
was  measured  and  calculated.  This  factor  of  2 by  itself,  of  course, 
does  not  account  for  the  difference  between  the  calculated  results 
using  our  program  and  the  experimental  data  furnished  by  ETEC. 
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TABLE  5 

RELATIONSHIPS  OF  EDGE  SCAN  VERSUS 
% CURRENT  DIAMETERS 


Edge 

^meas 

% Current  Radius 

to 

Scan 

^total 

Measured  Radius 

x%  to 
l-x% 

^80/r 

' meas 

'^90/r 

meas 

'^99/r 

meas 

1-99 

0.932 

0.77 

0.93 

1.31 

5-95 

0.740 

1.10 

1.31 

1.85 

10-90 

0.563 

1.40 

1.67 

2.36 

15-85 

0.413 

1.74 

2.08 

2.95 

20-80 

0.302 

2.11 

2.53 

3.58 

Therefore,  the  search  continued  for  discrepancies  between 
the  calculated  and  experimental  results.  Another  discrepancy  was 
found  in  the  mode  of  operation  of  the  two  lens  system.  Our  two 
lens  analytical  programs  are  set  up  for  the  special  case  of  a 
parallel  beam  between  lenses.  The  reason  for  this  is  that  it  is 
well-known  that  the  aberration  coefficients  for  a given  lens  are 
polynomial  functions  of  the  magnification.  Therefore,  the  lowest 
aberration  coefficients  occur  when  the  lenses  are  operated  with 
either  the  image  or  the  object  at  infinity.  A design  optimized 
for  this  condition  therefore  should  result  in  more  current  per 
unit  spot.  In  talking  with  Veneklassen,  he  indicated  that  his 
first  lens  was  optimized  for  the  thermionic  emission  gun;  however 
the  final  lens  was  the  objective  lens  from  an  ETEC  autoscan  SEM 
and  not  necessarily  optimized  for  the  field  emission  cathode. 
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His  analysis  therefore  indicated  that  running  the  beam  with  either 
a slightly  diverging  or  converging  beam  between  lenses  gave  more 
current  per  unit  spot.  In  order  to  study  this  effect,  some  of 
our  existing  programs  were  modified  to  allow  calculation  for  a 
nonparallel  beam  between  lenses.  At  present,  all  the  CR&D  optical 
programs  are  for  electrostatic  lenses  and  therefore  the  exact  case 
of  a first  electrostatic  lens  and  second  magnetic  lens  of  the 
Veneklassen  probe  could  not  be  treated.  We  were  able,  however, 
to  choose  two  electrostatic  lenses  that  came  close  to  the  properties 
of  the  electrostatic  and  magnetic  lens  of  the  Veneklassen  system 
for  analytical  purposes.  In  this  manner  we  could  study  the  effects 
of  operating  with  the  beam  parallel  or  nonparallel  between  the 
two  lenses.  The  results  of  beam  current  versus  spot  size  calcula- 
tions for  this  all  electrostatic  lens  system  are  shown  in  Figure  62. 
Shown  on  this  figure  are  cases  where  the  beam  is  parallel  between 
lenses  and  slightly  diverging  and  slightly  converging  for  the  case 
when  the  final  lens  is  not  optimized  with  the  first  lens.  Also 
shown  on  this  curve  is  parallel  beam  case  where  the  final  lens 
was  optimized  using  the  techniques  developed  in  previous  ion  column 
design  work.  From  the  figure,  it  can  be  seen  that  for  a nonoptimized 
lens  system  a slightly  diverging  or  converging  lens  yields  more 
current  per  unit  spot  than  the  parallel  beam  case.  This  is  in 
support  of  Veneklassen ' s claim.  However,  the  curves  also  show 
that  if  the  final  lens  is  optimized  for  parallel  beam  operation 
with  the  first  lens  then  a significant  increase  in  current  per 
unit  spot  can  be  achieved  with  this  mode  of  operation.  This  pro- 
cedure bears  out  our  previous  concept  of  the  proper  mode  of  opera- 
tion for  a two  lens  system.  These  calculations  show  that  another 


factor  of  two  in  beam  current  per  unit  spot  size  in  the  tenth  micron 
spot  size  region  could  be  achieved  by  running  the  nonoptimized 
lens  system  with  the  beam  slightly  diverging  or  converging  between 
the  two  lenses. 

In  further  discussions  with  Lee  Veneklassen  it  was  found  that  i 

i 

there  was  a difference  in  the  definition  of  the  chromatic  aberration  ! 

coefficient  between  his  usage  and  our  program.  Veneklassen  defines  | 

the  chromatic  aberration  coefficient;  | 

■■c  ' » (431 

whereas  the  definition  of  chromatic  aberration  coefficient  used 
in  our  programs  is; 


i 


V “ 


(44) 


As  can  be  seen,  the  chromatic  aberration  coefficient  is  a factor 
of  two  larger  by  our  definition.  What  this  means  is  that  in 
using  the  chromatic  aberration  coefficient  supplied  to  us  by 
Veneklassen  it  had  to  be  divided  by  two  in  order  to  fit  the  defi- 
nition for  our  program.  This  could  prove  significant  since  as 
can  be  seen  from  the  calculated  beam  current  versus  spot  size 
curves,  the  beam  is  operating  in  the  chromatic  limit  in  the  tenth 
micron  region. 


All  of  these  correction  factors  were  incorporated  in  our  op- 
tics analysis  program  and  the  beam  current  versus  spot  size  for 
the  Veneklassen  optical  data  was  recalculated.  The  results  are 
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BEAM  CURRENT.  Amps 


shown  in  Figure  63.  The  first  curve  shows  the  resulting  correc- 
tion due  to  the  factor  of  two  in  the  chromatic  aberration  coeffi- 
cient. The  next  higher  curve  shows  this  factor  for  along  with 
the  scaling  for  the  15-85%  rise  time  definition  of  spot  size.  Fi- 
nally, the  maximum  current  curve  shows  11  of  the  corrections  in- 
corporated; that  is,  the  factor  of  two  in  chromatic  aberration  co- 
efficient, correction  for  the  15-85%  rise  time,  and  finally  the 
scaling  of  the  factor  of  two  for  a nonparallel  beam  operation. 

As  can  be  seen  from  this  curve,  approximately  one  microamp  of  cur- 
rent is  obtained  in  a tenth  micron  spot.  This  is  in  excellent 
agreement  with  the  experimental  data  previously  shown  in  Figure  59 
for  the  Veneklassen  probe. 

SUMMARY 


During  this  reporting  period  negotiations  began  with  ETEC 
Corporation  for  a high  current,  thermionic-aided  field  emission 
microprobe.  Analytical  studies  were  made  in  order  to  better  under- 
stand the  column  specifications  and  performance  claims.  The  first 
analytical  attempts  were  grossly  in  error  with  ETEC's  published 
experimental  results.  Investigations  of  these  discrepancies  showed 
that  they  were  the  result  of  differences  in  definitions  of  spot 
size,  chromatic  aberration  coefficient,  and  a different  operating 
mode  (i.e.,  noncollimating  between  lenses).  Once  these  differences 
were  recognized  and  incorporated  in  the  analysis,  the  experimental 
and  calculated  beam  current  versus  spot  size  curves  were  in  excel- 
lent agreement.  As  a result  of  these  calculations  we  have  extended 
our  confidence  in  the  electron  optical  analytical  techniques. 


115 


1 

Also,  we  have  a clear  picture  of  the  specifications,  expected  per- 
formance, and  potential  for  further  improvement  in  the  ETEC  column. 

Based  on  the  optics  specifications,  and  the  confidence  that  these 
will  be  met  or  exceeded  through  the  use  of  the  analytical  tech- 
niques, an  order  was  placed  with  ETEC  for  the  high  current  TFE 
microprobe.  < 


APPENDIX 


H0TSP0T3 

DESCRIPTION  AND  USER’S  GUIDE 


INTRODUCTION 


HOTSPOT3  is  a finite  difference  type  computer  program  designed 
to  permit  the  study  of  heat  production  and  flow  in  any  solid  mate- 
rial in  which  the  assumption  of  radial  symmetry  for  both  the  mate- 
rial and  temperature  distribution  is  valid.  It  is  physically  a 
three-dimensional  solution  which  has  been  reduced  to  two  dimensions 
mathematically  by  assumption  of  radial  symmetry.  In  practice,  a 
solid  object  is  'constructed*  of  various  materials  whose  relevant 
properties  are  described  in  functions.  The  initial  conditions  are 
described  by  the  user  along  with  the  manner  in  which  heat  will  be 
added.  The  program  then  prints  on  an  output  file  an  array  showing 
the  temperature  distribution  everywhere  in  the  solid  object:  that 
is,  the  temperature  of  each  cell  at  any  future  time  chosen  by  the 
user . 

The  program  described  here  is  written  in  the  BASIC  language 
of  the  General  Electric  ISBD  tine  sharing  (Mark  3)  system. 


Program  Description 


1.  Mesh  Dimensions 

Since  the  problem  is  solved  in  cylindrical  coordinates,  the 
mesh  is  shaped  accordingly.  Thus,  the  volume  of  cells  at  small 
values  of  R is  smaller  than  that  at  large  R.  The  user  must  decide 
how  many  cells  are  needed  for  the  problem.  There  are  R9  cells 
in  the  radial  direction  and  Z9  in  the  axial  direction.  These  values 
are  set  in  lines  1200  and  1300  and  the  dimension  statement  of  line 
900  may  need  changing.  The  cell  lengths  are  R1  and  Z1  in  units 
of  (cm)  and  are  set  in  lines  1600  and  1800  (cm) . This  version 
of  HOTSPOT  is  intended  to  model  a substrate  with  up  to  two  uniform 
surface  layers.  The  top  layer  consists  of  cell  row  numbers  1 
through  (Z2-1) . The  second  layer  begins  with  row  Z2  and  ends  with 
row  (Z4-1) . The  substrate  rows  are  Z4  through  Z9.  Z2  and  Z4  are 
set  in  lines  1700  and  1750.  It  is  generally  wise  to  first  try 
a problem  with  rather  a coarse  mesh  with  small  R9  and  Z9  until 
the  sort  of  flow  pattern  and  time  scale  is  understood.  The  cost 
increases  rapidly  with  these  parameters. 

2.  Image  Matrix 

The  image  matrix  I(R,Z)  is  an  array  corresponding  to  the  cells 
of  the  mesh.  The  elements  of  I(R,Z)  are  integers  denoting  the 
material  each  cell  is  made  of.  The  matrix  is  formed  beginning 
on  line  3000. 
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3.  Power  Input 


The  total  input  power,  P,  is  set  in  line  1100.  P is  given 
as  a function  of  EO,  the  bezun  energy  in  kilovolts,  with  a one  micro- 
ampere beam  current  assumed.  This  may  be  changed  at  the  discretion 
of  the  user.  The  only  remaining  parameter  left  to  the  user  is 
the  beam  radius  (Gaussian  width  at  the  surface) , Al,  set  in  line 
1030.  Even  though  the  beam  power  function  is  defined  in  absolute 
terms,  the  result  is  renormalized  so  that  the  sum  of  the  blocks 
of  power  supplied  to  all  the  cells  is  equal  to  the  specified  input 
power.  The  power  density  function  FNP(R,Z)  is  formed  beginning 
at  line  12900. 

4.  Time  Scale 

Each  time  increment  must  be  chosen  small  enough  so  that  only 
a relatively  small  change  occurs  in  the  temperature  of  the  cells 
or  the  program  will  not  give  correct  results.  This  minimum  time 
interval  T1  is  set  in  line  2100  and  its  optimal  value  depends  on 
the  physical  size  of  the  volume  being  heated  and  the  thermal  con- 
ductivity of  the  materials.  The  way  the  program  operates  is  to 
form  an  increment  cell  matrix  D(R,Z)  whose  elements  consist  of 
the  net  temperature  change  for  each  cell  that  occurs  during  a time 
T1  based  on  its  temperature  and  that  of  its  neighbors.  This  matrix 
is  then  added  to  the  temperature  matrix  U(R,Z).  This  process  con- 
tinues until  a time  T9  (set  in  line  2000)  is  reached  when  the  ma- 
trix U(R,Z)  is  stored  on  an  output  file  with  suitable  headings. 


i 

I 

I i 

The  program  will  then  ask  for  a new  value  of  'T  MAX'  (T9)  and  if 
given  it,  will  proceed  as  before.  If  it  is  desired  to  stop  the 
calculation  but  retain  the  option  of  continuing  later,  respond 
with  a carriage  return  and  all  of  the  data  will  be  saved  in  the 

J 

file  'HEAT2'.  To  continue,  copy  'HEAT2'  into  the  input  file  'HEATl' 
and  RUN.  If  for  any  reason  you  do  not  wish  to  have  HEAT2  altered 
at  the  conclusion  of  a run,  type  'S'  and  the  program  will  stop. 

At  each  stop,  the  results  at  that  time  (T9)  will  be  printed  in  | 

the  output  file  'HEATOUT'.  This  file  is  always  appended  (never 
destroyed)  by  the  program.  When  run,  the  program  first  asks  'NEW 
OR  OLD?'.  If  the  answer  is  'NEW',  it  will  start  with  a clean  slate 
from  time  T=0.  Otherwise  it  will  read  'HEATl'  and  use  those  data 


for  starting  conditions. 

5.  Boundary  Conditions 

As  the  program  stands  at  present,  all  of  the  outermost  side 
and  bottom  cells  are  maintained  at  the  substrate  temperature  (set 
in  line  5200)  which  is  also  the  initial  temperature  when  a 'NEW' 
case  is  run.  These  are  all  cells  with  R=R9  or  Z=Z9.  At  present, 
the  Z=1  cells  are  not  influenced  by  the  outside  environment  (except, 
of  course,  by  the  input  power  function).  Provision  is  made,  how- 
ever, to  allow  heat  to  flow  in  or  out  of  the  top  surface  (e.g., 
radiation)  in  line  8400. 
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6.  Volume  Specific  Heat 


The  volume  specific  heat  function  FNC(T,I)  allows  a tempera- 
ture dependent  specific  heat.  I is  the  material  index  for  a par- 
ticular material  specified  originally  in  the  image  matrix.  Each 
cell  of  the  mesh  is  assigned  a number  which  specifies  the  material 
the  cell  volume  is  made  of  (see  Image  Matrix) . 

7.  Thermal  Conductivity 

This  function  FNK(T,I)  is  defined  following  line  16500.  Again, 
provision  is  made  for  temperature  dependence. 

8.  Density 

The  material  density  is  used  in  defining  electron  range  in 
materials  when  one  is  constructing  the  power  input  function.  This 
function  FNR(I)  is  a table  of  values  whose  definition  begins  in 
line  19600. 

I 

9.  Material  Label 

I 

I 

This  is  a string  function  FNL$(I)  that  returns  the  name  of 
the  material  designated  by  the  index  number  I.  Its  use  is  for 
* automatic  printouts  of  materials  used  in  each  calculation. 

I 


I 
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10.  Output 


The  output  to  the  file  HEATOUT  begins  with  the  writing  of 
the  temperature  matrix  U(R,Z)  as  it  would  appear  after  a time  T1 
in  the  absence  of  heat  flow.  This  is  useful  to  help  visualize 
the  effect  of  the  pattern  of  input  power.  Just  preceding  this, 
and  in  every  printing  of  U(R,Z),  is  a set  of  headings  written  by 
lines  9300-10000  and  alterable  at  the  discretion  of  the  user. 

The  first  of  these  lines  prints  the  comment  string  C$  defined  in 
line  800.  The  others  help  identify  the  case  being  run. 

In  addition  to  the  output  file,  there  is  output  to  the  ter- 
minal which  is  useful  for  observing  progress  as  time  increases. 
This  printout  is  due  to  lines  10900-11500  and  at  present  is  part 
of  the  first  row  and  column  of  U(R,Z);  enough  to  guide  selection 
of  time  intervals  without  taking  excessive  time  at  a slow  termi- 
nal. These  lines  are  also  alterable  according  to  the  wishes  of 
the  user. 


HOTSPOT3  USER'S  GUIDE 


1.  Preliminaries 

The  first  step  is  to  set  the  values  of  the  various  parameters 
in  the  progreun  to  match  the  case  being  run.  For  most  parameters 
this  is  simple,  but  a few  are  tricky.  The  most  difficult  is  Tl, 
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the  time  of  one  iteration.  It  appears  from  our  studies  that  one 
should  use  the  largest  possible  value  of  T1  for  which  the  process 
is  stable.  Use  of  a smaller  value  does  not  substantially  change 
the  result.  A good  first  try  would  be  about  1%  of  the  cell  area 
divided  by  the  diffusivity  although  this  can  vary  widely  with  the 
geometry  of  the  problem.  It  is  our  practice  to  first  estimate 
Tl  in  the  manner  described,  setting  T9  to  ten  times  this  value, 
and  then  run  the  program.  If  the  process  is  unstable  as  evidenced 
by  wild  oscillations  in  the  cell  temperatures,  reduce  Tl  by  half. 
Otherwise  double  it,  repeating  this  procedure  until  the  largest 
stable  value  is  found.  Decreasing  the  cell  size  R1  and  Z1  gener- 
ally requires  decreasing  Tl.  It  is  recommended  that  the  cell  size 
be  made  as  large  as  reasonable  at  least  for  preliminary  runs. 

It  is  necessary,  however,  that  each  separate  material  layer  have 
at  least  two  cell  rows.  The  following  is  an  alphabetic  listing 
of  all  of  the  simple  parameters  which  would  normally  be  changed 
by  the  user  (with  line  number) ; 

A1  (1030)  Beam  Radius  (cm) 

EO  (1010)  Beam  Energy  (kilovolts) 

P (1100)  Beam  Power  Input  (watts) 

R1  (1600)  Radial  Cell  Dimension  (cm) 

R9  (1200)  Number  of  Cells  in  the  Radial  Direction 
R9  (900)  Check  DIM  of  V,M,U,P,D,I 

Tl  (2100)  Time  of  One  Iteration  (sec) 

T9  (2000)  Time  of  First  Program  Stop  (sec) 
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Zl  (1800)  Axial  Cell  Dimension  (cm) 

Z2  (1700)  Index  I of  First  Row  of  Second  Layer 
Z4  (1750)  Index  # of  First  Row  of  Substrate 
Z9  (1300)  Number  of  Cells  in  the  Axial  Direction 
Z9  (900)  Check  DIM  of  V,U,P,D,I 

Next,  the  image  matrix  I(R,Z)  must  be  constructed  to  agree 
with  the  configuration  being  modelled.  At  present,  the  program 
has  a knowledge  of  7 materials.  To  discover  what  number  corre- 
sponds to  a given  material,  read  the  function  subprogram  for  the 
label  function,  FNL$(I),  beginning  on  line  18500.  If  the  desired 
material  is  not  listed,  it  will  be  necessary  to  add  it  as  described 
below.  If  the  program  is  to  be  used  with  only  one  layer  over  the 
substrate,  give  the  middle  layer  the  number  of  the  substrate  or 
delete  line  3720. 


Line  3300:  defines  top  layer  material  I 
Line  3720  defines  middle  layer  material  # 
Line  3600  defines  substrate  material  I 


The  string  C$  defined  in  line  800  is  an  identifying  comment 
which  will  be  printed  with  the  output. 


When  it  becomes  necessary  to  add  a new  material,  change  a 
material,  or  vary  the  data  for  an  existing  material.  The  material 
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related  function  subprograms  must  be  changed.  The  arguments  T 
and  I are  temperature  and  material  number.  The  four  functions 
are: 


FNC(T,I)  (14600)  Volume  specific  heat  ( joules/cc-deg  K) 
FNK(T,1)  (16600)  Thermal  conductivity  (watts/cm-deg  K) 

FNL$(I)  (18500)  Name  of  material 
FNR(I)  (19700)  Density  (gm/cc) 

2.  Program  Output 

As  discussed  earlier,  the  program  writes  results  on  the  file 
HEATOUT  as  well  as  the  terminal.  An  enormous  quantity  of  informa- 
tion is  calculated  but  usually  only  a limited  amount  is  preserved. 
The  full  array  of  temperatures  is  calculated  at  each  interval  of 
time,  Tl.  This  is  the  principal  output  of  the  program  and  gives 
the  temperature  of  each  cell  at  that  time.  The  two  rows  of  tem- 
peratures printed  below  the  array  are  the  temperatures  at  the  in- 
terface of  each  of  the  layers.  In  alloy  studies,  these  are  the 
most  important  results.  It  is  up  to  the  user  to  decide  at  which 
times  the  temperature  array  should  be  recorded.  The  first  record 
is  made  at  time  T9  set  by  the  user.  After  that,  the  progr2un  will 
ask  for  the  next  time  that  a temperature  record  should  be  made 
(NEW  T MAX?).  Even  if  only  a 'final'  temperature  is  desired,  it 
is  useful  to  plot  temperature  as  a function  of  time  to  verify  that 
a near  final  temperature  has  been  reached.  Normally,  frequent 


stops  should  be  made.  It  is  the  last  opportunity  to  obtain  these 
data. 


3.  Operation 

The  effective  use  of  the  data  storage  files  HEAT2  and  HEATl 
requires  some  thought.  At  every  stopping  point,  all  current  data 
are  written  into  the  file  HEAT2  unless  the  character  'S’  is  used 
as  a response  to  a request  for  the  next  time  (NEW  T MAX?) . The 
purpose  of  this  feature  is  to  permit  recovery  of  data  in  the  event 
that,  alter  a long  calculation,  the  time  interval  Tl  should  need 
to  be  shortened  (as  occasionally  happens) . If  this  should  happen, 
the  temperature  data  at  that  time  would  be  useless  and  typing  'S' 
would  prevent  the  destruction  of  the  previously  stored  data.  The 
program  would  then  stop  permitting  a change  in  the  value  of  Tl 
and  could  then  be  run  again.  At  this  running,  the  word  'OLD*  should 
be  typed  in  response  to  'NEW  OR  OLD?'.  The  last  leg  of  the  calcu- 
lation will  be  repeated  using  the  new  value  of  Tl.  In  the  event 
of  a normal  termination,  the  final  data  will  be  stored  in  HEAT2 
where  it  may  be  used  to  extend  the  calculation  or  it  may  be  stored 
under  a different  name  for  possible  future  use.  Remember,  as  dis- 
cussed earlier,  the  contents  of  HEAT2  must  be  moved  to  HEATl  to 
be  used.  HEAT2  is  always  an  output  file  and  HEATl  is  always  an 
input  file. 
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The  following  is  a listing  of  the  program. 

10  ' P.UhAV/H«fSP0T3 
20  ^ 

30  ' SPECIAL  VEhSlON  OF  HOTSPOT.  FIXED  HI, 7.1,  Tk<0  UNiFOhM  FILMS 

40  ^ ON  substhate.  complete  electron  beam  algorithm  built  in. 

DO  ' FOR  information  CALL 

00  ^ PETEk  V GRAY 

70  ' d*235  P369  Ok  (SIB)  3R5  8369 

dO  * 

100  CALL  k.AHN(n) 

300  FILES  FEATOUr 

400  APEENi;#! 

bno  MARGIN#! , 124 

oOO  FILES  HEATUHEAT2 

/OO  REST0RE#2 

800  C$="COMMENT  LINE" 

900  DIM  V(!5,l5),Ki5,2),U(l5,!b),P(l5,i5),D(l5,!5),I(l5,l5) 

910  DIM  L(20) 

1010  E0=4  'BEAM  KILOVOLTS 

1020  Fl=l  'FLAG  FOR  FNP 

1030  AI«500E-8  'BEAM  RADIUS 

lino  p=E0*lE-3 

1200  R9»I5 

1300  Z9=I5 

1400  R9*R9-1 

1500  Z9=Z9-1 
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1600  kI*200E-8 


1700  Z2=4 

1750  Z4»6 

1600  Zl=200E-8 

1900  T=0 

2000  T9=2.bE-ll 

2100  T1=2E-12 

2200  FkINT'‘NEW  OH  OLD"  * 

2300  INkUTAS 

2400  IFAS  = ''NEk"THEN3O00 

2600  MATkEAD#2,U,P,D,I 

2700  RbAD»2,T 

2750  RtAU*2,C$ 

2800  GOTO  11 700 
2900  ' 

3000  'FORM  IMAGE  MATRIX 

3100  F0RR=1T0R9+1 

3200  F0HZ=1T0Z2-1 

3300  I(R,Z)=6  'TOR  LAYER 

3400  NEXTZ 

3500  F0RZ=Z2T029+ I 

3600  I(R,Z)=1  'SUBSTRATE 

3700  NEXTZ 

3710  F0RZ=Z2T024-1 

3720  I(R,Z)«7  'MIDDLE  LAYER 
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3730  NEXTZ 
3800  NEXTI< 

3900  ' 

4000  ^FORf*  POrtEh  1-ENSITV  MATRIX 
4100  Pl»0 
4200  F0RR»IT()R9 
4300  F()HZ»IT0Z9 

44  00  KH,Z)«FNF(<k-.5)*RI , (Z-.5)*ZI ) 

4bO0  PI*PI+P(R,Z)*3. I 4)6*ZI*Rl*Rt*(2*R-I ) 

4600  NEXTZ 
4700  NEXTR 
4800  MATP»(P/F1)*I 
4850  PRINT''PI«*'PI 
4H60  PRINTC$ 

4864  PRINT''A*'(Z2-1  )*ZI*1E8i"ANGSTR0M  ••  f FNLS  ( I ( 1,1  ) ) * *'  FILM  ON  "I 
4868  IF  I ( I ,Z2)*I ( 1 ,Z9+1 ) THEN4RR0 
4872  print 

4676  PRlNT“A«(Z4-Z2)*ZI*IE8f  "ANOSTROM  "IFNLSI  K J ,Z2  )) »"  FILM  ON  •'» 
4«80  8RINTFNL$(I( J,Z9+I >)l"  SUBSTRATE." 

4684  PRINT"TIME«"T 

4PFR  PKINT'E0="E0'’KV,  P()rtER»''P''RATrS" 

4R92  PRINT''REAM  RADIUS="Al  "CM" 

4896  PRINT" VOL  SP  HEAT*"FNC(T,I ( I , D)  lFNC(T, I ( I ,Z2 ) ) I 
489/  PRINTFNC(T,I(I,Z9+1 )) 

4900  PRINT"TH  COND»"FNK(T, I ( 1,1 ))  IFNK (T, I (4  ,Z2 ) ) IFNK (T, I ( 1 ,Z9+1 )) 
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4904  FhINT''DRa“Wl  l''DZ=''Z1 


4903  f-hINT 
boon  'THfcPMAL  BIAS 
5100  .V.ATU*CON 
5?<^0  WATU=(300)*U 
5300  ' 

5400  'FIRST  lOvVEH  TO  U(H,7.) 

5500  F0Wk=IT0R9 
5600  F0hZ=IT0Z9 

oVOO  U(R,Z)=U(R,Z)+TI*P(H,Z)/FNC(U(R,Z),I (R,Z)) 

5PO0  NEXTZ 
d900  NEXTR 

6000  FklNT#1  ,*'HhAT  ADDED  DURING  FIRST'*  »T I l" SECONDS'* 

6100  FhINT'*HEAT  ADDED  DURING**lTl  "SEC** 

6200  GOT09300 
6300  T=TI 
6400  ' 

6500  'FOR?.,  INCREMENT  MATRIX  D(R,Z)  INCLUDE  FLON  AND  POl'.ER 
6600  f0RR=lT0R9 
6700  FORZ=ITOZ9 
6b00  U=U(H,Z) 

6900  W*I (k,Z) 
iOnO  K*FNK(U,M) 

/mo  C«FNC(U,M) 

7200  A2«K/C 
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i3C0  l.=Tl*k(H,Z)/C 
/^ro  I I THEN  7700 

/boo  i'=D+/s2*TI*((P-l  )*U(K-I  ,Z)-(2*k-l  ) *U+R*U ( h+ 1 ,Z  ) )/ ( < W-.6)*RI  *R  I ) 
Zero  OOTO'/C' no 
t(ro  ■'psl  case 

liiro  b=|'+/!2*Tl*2<f{U(2,Z)-’J)/(Rl*RI  ) 

/<,0()  IEZ=ITHbNH20n 
/vin  lrUn,Z)<>UH,Z+l  ) THbmopn 

/>?r'  lEi  (i;.z)<>l  (k,z.-i  )THENao40 

wroo  DsO+AZ^TUCUCh.Z-l  ) -?>^U+U  ( R , Z+ 1 ) )/(ZI*Zl  ) 
bOlO  uf-'TOf.-bop 

bC20  U=U+A2*1  U(o(H,Z-l  ) -2*U+V  ( K , Z+ I ))/(ZI*Z.I  ) 

.5030  uC.TO'  boO 

oOAD  L=i>A2*Tl*(  V(P,Z-I  )-2*U+b{R,Z+l  ))/(ZI*/lI  ) 
a 100  uOrObbno 

•5200  '^.=  \ CASE—CC^Nbl  Ct.h  KAOIATK'N  LOSS 
ojOO  u=r+/.2*TI*(U(R,2)-U)/(Zl*Zl  ) 
oA"-'-'''  'L'=h+PADIAT  ION  L('SS 
'■5^00  i (W,7)=[; 

Mcro  .ihXTZ 
r,  <n  UhXlM 
'{•on  ' 

0900  1=T+T1 
yO'^'O  MATU=u+D 
yrO'3  'rOh/:.  .MIRROR  f.'ATKiX 

i 
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9CI0  H()RR=IT0f<9 


9015  KI=HKK(U(R,Z2-I ) ,I (K,Z2-I )) 

9020  K2=FNK(U(H,Z2) ,I (R,Z2)) 

9025  UI=U(k,Z2-l ) 

9030  U2=U{R,Z2) 

9035  ,V,(R,  I ) = (U1*KI+U2*K2)/(.KI+K2) 

9040  V(n,Z2-I )=2*M(H, 1 )-U2 

9042  V(R,Z2)=2*M(k, 1 )-UI 

9055  K I =FI  'K  ( U ( H , Z4- 1 ) , I ( H , Z4- 1 ) ) 

9060  K2-hNK(U(R,Z4) ,I (R,Z4)) 

90o5  U1=U(R,Z4-I ) 

9070  U2=U(k,Z4) 

9075  rt ( K , 2 ) = ( U I I +U2*K2 ) / ( K I +K2 ) 

9080  V(rt,Z4-l)=r2*J/(k,2)-U2 
9085  V(rt,Z4)=2*fr(H,2)-Ul 
9090  NEXTh 
9095  IhT=0THEN6300 
9100  ' 

9200  IET<T9-T1/2THEN6500 
9300  HklNT#l,C$ 

9400  FklNT#l  ,*’A“(Z2-1  )*ZI*IE8»''ANGSTR0M  " lENLS  ( I ( I , I ) ) i ••  FILM  ON  "* 
9420  IF  I( 1 ,Z2)=I ( 1 ,Z9+I ) THEN9460 
9440  FRINT#I , 

9450  FRINT#I  ,"A"{Z4-Z2)*ZI*IE8l''ANGSTkOM  •' IFNI.S  ( I ( I , Z2  ) ) I 

9451  FRINT#I  FILM  ON  "I 


^ \ •4*** 


. > ' u.»- 


y4oO  K.lNfrl  ,FffLS<i  (I  ,4V+n)  t"  HUMSTf?/Tc 
y^po  I Ml 1 , "TIi  'iF  = " T 


II 


ycoo  f-rlin^  I ,“=r'=**HO"r.V,  A'rrS'* 

HhlKT-' 1 ,-HLAA  h Ai 'TuS="  A I ••(:.'.'•• 

y/PT  r-FlNln  I VOL  SP  !h  AT="h  HC  (T,  I ( I , I ))  IFNC  C T,  I ( I , Z2 ))  I 
y <01  f kINTv^I  ,P.MC(  r,I  ( 1 ,Z9+1  )) 

yPOO  PRIK1>1  ,"TH  C()ND=“FNK(T,  I (1,1  ))  tFNraT,  I ( I ,72))  J 
ycO)  FhiiiT^^l  ,FNK(T,I  ( 1 ,Z9+n) 
y<rPO  FKi.NTf  I ,*'Dh='*RI 
IC'OOO 

I 0100  F()«Z=IT0Z9 
10200  l-()RF=IT()riV 

K^j'^O  PRI:n^*l  ,TAJ(b*(P-l  ))  »STPS(U(K,Z))  » 

10400  NhXTP 
I05'''0  FhII!T#l, 

10600  MtXTZ 
10(00  IkI'T#!, 

10710  H()h/  = 1T02 
10720  F()W1;=1T()h9 

IOVjO  HRINT((1  ,TA.^(6*(k-I  ))  lSTP$  (/' ( h,  Z ) ) » 

10/40  MEXTn 
lOVbO  PKiNlill, 

10)60  NFXTZ 
10 //O  ^KINT#I , 

lOHOO  Fr<IM</|  , 
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10900  F()PP=IT()W9 
I 1 ono  Hm1NTU(R, 1 ) I 
I I I oo  NhXTW 

I 1200  FkIhT 

II  300  F()R^.=  2T0Z9 

1 Moo  ^RINTU(  1 ,Z) 
llDOO  NbX'U 

I l:;10  l-HIKT“b()UNUKY  ThWFHWATUkES="M  ( I ,1  ) tM(l  ,2) 
lloOO  IrT=0THEN9005 

II  Vf'O  FrtIM"T=''Tl",  NLt'.  T .VAX"  I 
IlftOO  rATlNPUT  'fS 

12fiO:)  Ib'IS(  1 )=“S"THhN12/no 

12300  SCKATCH4i^3 

12t)O0  .’A1‘  bRITE  «-‘3,U  JKlUlI  » 

126O0  ,u(’IThn^3,l 
1 2CD')  v.Hl'lt  .'t 
12  660  IFNUA=0:Hbil  1 2 /on 
126/0  1V  = VAL(TS( 1 )) 

I2c,'v0  00106600 
12/00  bTOF 
120  00  ' 

1 2900  'hOKtU  DENSITY  FUNCTION  (UHNORMAl.IZtU) 
IjOOO  Lit;FFIiF(X,Y)  't.ATTS/CC 
13010  'VthTlCAL  eUflCTlON 
13013  IF  FI=0TliEi'il3090 


13020  'FIRST  TIME  THROUGH 
15030  F0RI=0T()I3 
13040  hEAD  L<I) 

I30b0  DATA  .75, 1 .05, 1 .27, I .4, 1 .4 j ,31 J .08» .82, .55, .32, 

13051  DATA  . 18,. 08,. 02,0 
13060  NEXTI 

13070  Z3»(Z2-1)*Z1  'METAL  THICKNESS 

13030  Z0=Z3*(FNR(I(1,1 ))/FNH(I(l,Z2))-l)  'DISPLACEMENT  OF  ORIGIN 
13(>82  Z5*Z0+Z3+(Z4-Z2)*Z1 
1 3034  Z6-Z5* ( FNW ( I ( 1 , Z2 ) ) /FNR ( I ( 1 , Z9 ) ) - 1 ) 

13035  F1*0  'RESET  FLAG 
IJOVO  'E7ERY  TIME  THROUGH 

13095  A=A1+Y*TAN(3. 14)59*30/180)  'mIDTH  AT  1 /E  POINT 

131O0  IFY<Z3THEN1 3120 

13110  Y=Y+ZO  'FORM  EFFECTIVE  Y 

13112  IFY-Z0<(Z4-1 )*Z1THEN13120  'MIDDLE  LAYER  CASE 
13114  Y=Y+2o  'SUBSTRATE  DISPLACEMENT 

13120  kO=EO''I  .75*4E-6/FNR(I  (R,Z))  'GRUN  RANGE 
13125  Y*Y/R0 
13130  L=0 

13140  IFY>=1  .3THcM3180 
13150  Y1»1NT< 10*Y) 

13160  Y2»10*Y-Y1 

131/0  L-L(Y1  )-»-Y2*(L(Yl  + l )-L(Yl  )) 

13180  V«L*E0*lE-3/R0  'VERTICAL  PART 
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'liOnlZONTAL  hUF!CTI()M 
1 i2  1 0 F-  = £:XP(  -X*X/  ( A*A  ) ) / ( A*A*  J .14)59) 
13220  rNP=H*V 
Ij2oO  PNtM. 

144ri')  * 

14500  'VOLDf'FE-:  S^'cCl^■IC  HtAT 

Uo'^O  l.hFPNC(T,f.')  'JOULES/CC-rtCJ  K 

14/r'O  IHM=ITHENIb  I'Xi'bl 

MHOO  JF,V.=2ThENl3400‘'GH 

I4VO0  lhP.=3TI-EMIt> /OO'AL 

1500,')  It-M=4THbN'loOnO‘'AU 

15010  Ihf.'.=5THENIo2IO‘''lE 

15020  IFM=6THEN1o225  'SELENIUM 

I 50J0  1 E .'.i=VTH£N  1 0240  ' 5102 

15100  'SILICON 

I 5200  i-NC=I.63 

15300  001016300 

15400  ^GEkMANIU/.', 

15500  ENC=I.65 
156(X)  0010I0300 
15/00  'ALUMINUM 
150O0  i-NC=2.43 
lt;S?00  00  re  I 6300 
I 6 000  'oOLD 
16100  ENC=2.5 
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IA2'^0 

GOT0 16300 

16210 

'TELLURIUM 

16216 

hNC=2.39 

1 622) 

GOT(>  16300 

1 6225 

'SELENIUM 

16230 

ENC=.94l 

lo235 

GOTO  1 6300 

16240 

'SI02 

16244 

RNC*1 .63 

16250 

GOTO  16300 

lo300 

FNEND 

16400 

/ 

l65rx) 

'THERMAL  CONDUCTIVITY 

1 ^600 

CEFFNK(T,W)  'rtATIS/CM 

16/00 

Ihw  = lTHENl  ■/ lOO'SI 

1 ct'OO 

IHM=2THEN1 74O0'GE 

1 6900 

IrM=3THEN1 7700'AL 

1 /0f>0 

I rM=4THEN 1 dOOO' AU=3 . 0 

1 7010 

IHM=5TIEN’13210  'TE 

1 7020 

IHo=oTHEN13?25  'SE 

1 /030 

I EM= /THEN  16240  'SI  02 

1 / iro 

'SILICON 

1 72'^0 

HNK=.d3 

1 /3O0 

GOTO  18300 

1 y4fX) 

'GERMANIUM 
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••  * -r  * 

* • J»  ■ ^ -»-i  • * 

• V.  ' ■ 


1 FNK=.o^ 

I /ooo  o0T()l8300 
IV/^'O  ''ALu*' I.sU.’.l 
I /8'io  rNN=.4V" 
l/v'-O  U)1('I.m3i^O 
IrOf'.i)  ''uOLD 
U K')  ri'.K=3 
ItPaO  u0T('K3'^0 
lf2i)  ''TfcLLUFIuM 
U 2 I 5 ^ 'tK- • 06 

U)J(;1H3^0 
U-22‘J  'btLb.NlUV 
i;.2jO  F.'jN=.6n5 
lf-2o‘5  l()T('I83'^0 
I > 240  'M02 

\^-24-j  hllfv=.0lt3 

1'-2-jO  U)T('1<300 
If  Jl'O  FucI'L' 

U-4ro  ^'lAThwlAL  LAbtL 

ihhrivLS(]) 

lfc6M0  ():>i  1 GOTO  I H?on,  1^900,  K too,  19300,  19410,  19420«  19430 
It'KO  1-NL$=*'SILIG(H.'** 

1^8ro  U)T0l9bno 

1 i-9f'0  hifLS  = ''Gc9MAI.IUV,'* 


ivr'f.O  uO  TO  19b  00 


I ^100  FNL$***ALUM1NUM" 

I 92 no  L()TOI9bno 
l93no  FnLS=“oOLD“ 

I 9^ no  00X019500 
19410  h NL  S = “Th LI. J 1<  I U M “ 

I 94  It)  GOTO  1 9500 

19420  hNLS  = “SELbNIUM'' 

ly4  2D  oO 10 1 9500 

I94j0  hNL$='*SILICON  OlOXIOb'' 

19435  U)T0l95O0 
195O0  hNtNli 

l9oro  ^OcKSITY  FUNCTION 
19/00  LhrFNR(I)  'OM/CC 

I9BO0  ON  I GOTO  I 9900, 20 1 OO, 20300, 20500, 206 1 0, 20620. 20630 

19900  rNk=2.42  'Si 

20000  O0T0207O0 

20l''’'0  rNii=5«36  'Gh 

20200  U)T020/00 

2C300  FNh=2.699  'AL 

2(;4O0  U)T('20700 

2050'J  FNk=l9.32  'AU 

2O0OO  OOi()20VOO 

20610  FNH=6.24  'TE 

20615  G()T()20700 

20620  FNW«4.79  'SE 
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